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In this tutorial we review physical implementation of quantum computirig using a system of cold 
trapped ions. We discuss systematically all the aspects for making the implementation possible. 
Firstly, we go through the loading and confining of atòmic ions in the linear Paul trap, then we 
describe the coUective vibrational motion of trapped ions. Further, we discuss interactions of the 
ions with a làser beam. We treat the interactions in the travelling-wave and standing-wave config- 
uration for dipole and quadrupole transitions. We review different types of làser cooling techniques 
associated with trapped ions. We address Doppler cooling, sideband cooling in and beyond the 
Lamb-Dicke limit, sympathetic cooling and làser cooling using electromagnetically induced trans- 
parency. After that we discuss the problem of state detection using the electron shelving method. 
Then quantum gates are described. We introduce single-qubit rotations, two-qubit controUed-NOT 
and multi-qubit controlled-NOT gates. We also comment on more advanced multi-qubit lògic 
gates. We describe how quantum lògic networks may be used for the synthesis of arbitrary pure 
quantum states. Finally, we discuss the speed of quantum gates and we also give some numerical 
estimations for them. A discussion of dynamics on off-resonant transitions associated with a quali- 
tative estimation of the weak coupling regime and of the Lamb-Dicke regime is included in Appendix. 



PACS numbers: 03.65. Ud, 03.67.Lx, 32.80. Pj, 32.80.Ys 
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I. INTRODUCTION 

Although trappcd ions have found many applications in physics [Q, they caused a turning point in thc cvolution 
of quantum computing when the paper entitled Quantum computation with cold trapped ions was published by Cirac 
and ZoUer in 1995 This proposal launched also an avalanche of other physical realizations of quantum computing 
using different physical systems, from high finesse cavities to widely manufactured semiconductors Through the 
years we have learnt a lot, but also revealed many pcculiarities, about the physical realization of quantum computing 
which has led to many discussions concerning the conditions under which we could in principle implement quantum 
computing in certain quantum systems. 

Before we give the list of requirements for the physical implementation of quantum computing we will introduce 
the fundamental terminology to appear throughout this paper. We will foUow the definitions in Ref. Q. 

• A qubü is a quantum system in which the logical Boolean states and 1 are represented by a prescribed 
pair of normalized and mutually orthogonal quantum states labelled as |0) and These two states form a 
computational basis and any other (pure) state of the qubit can be written as a superposition 

|V) = «|0)+/5|1) (1.1) 

for some a and (3 such that jap + = 1. It can be shown that we may choose a — cosí? and [3 = e"^ sini?. 
A qubit is typically a microscopic system, such as an atom, a nuclear spin or a polarized photon, etc. In quantum 
òptics a two-level atom with a selected ground \g) and excited |e) state represents a qubit. Hence the notation 
\g) and |e) is used for the computational basis instead of |0) and |1). For instance, some qubits can serve for 
lògic operations or the storage of Information. Then we refer to lògic qubits. Some others can be used especially 
for sympathetic cooling of lògic qubits and we may call them cooling qubits. Some further qubits can be used 
as a quantum channel for transferring the Information between distinct lògic qubits and then we refer to them 
as to a quantum data bus. 

• A quantum register of size N refers to a collection of N qubits. 

• A quantum gate is a device which performs a fixed unitary operation on selected qubits in a fixed period of time. 

• A quantum network is a device consisting of quantum gates whose computational steps are synchronized in time. 

• A quantum computer (processor) can be viewed as a quantum network or a family of quantum networks. 

• A quantum computation (computing) is defined as a unitary evolution associated with a set of networks which 
takes a initial quantum state (input) into a final quantum state (output) and can be interpreted in terms of the 
theory of Information processing. 

For the moment we presume that the foUowing five requirements (termed Di Vincenzo 's checklist) should be met in 
order to realize quantum information processing on a quantum system |^ . Actually, there are two more requirements 
for the case of the transmission of qubits in space (flying qubits). However, it appears that all these requirements are 
necessary but not sufhcient for successful experimental realization of a quantum processor |^ . 

(1) The system must provide a well characterized qubit and the possibility to be scalable in order to create a 
quantum register. 

(2) We must be able to initialize a simple initial state of the quantum register. 

(3) Quantum gate operation times must be much shorter than decoherence times. The quantum gate operation 
time is the period required to perform a certain quantum gate on a single qubit or on a set of qubits. The 
decoherence time approximately corresponds to the duration of the transformation which turns a pure state of 
the qubit |V^) = a|0) + /3|1) into a mixture p = |a|2|0>(0| + |/3|2|1)(1|. 

(4) We need a set of quantum gates, to perform any unitary evolution operation that can be realized on the quantum 
system. It has been shown that any unitary evolution can be decomposed into a sequence of single qubit rotations 
and two-qubit controUed-NOT (CNOT) gates 0. 

(5) Thc rcsult of a quantum computational process must be eSiciently read out, i.e. the ability to measure distinct 
qubits is required. 
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Now we introduce briefly the physical system under consideration. Cold trapped ions is a quantum system of N 
atòmic ions confined in a linear trap. We assume an anisotropic and harmònic trapping potential. The ions are làser 
cooled to a very low temperature, beyond the Doppler cooling limit, reaching the recoil cooling limit Hence the 
term cold trapped ions. The ions form a linear crystal and oscillate in vibrational coUective motional modes around 
their equilibrium positions. In their internal structure, depending on the choice of atòmic species, we distinguish 
distinct atòmic levels. The ions are individually addressed with a làser or a set of làsers in the travelling-wave or 
standing-wave configuration. We can detect the internal state of ions using optical detection devices. Further, we 
address briefly the requirements for the physical implementation of quantum computing (mentioned above) using cold 
trapped ions. 

(1) The qubit is represented by a selected pair of internal atòmic states denoted as \g) and |e). This selection is 
discussed in detall in Sec. [V. The quantum register is realized by N ions forming the ion string in the linear 
trap, namely the linear Paul trap, which is reviewed in Sec.O. A selected collective vibrational motional mode 



(normal mode) is used as the quantum data bus. The vibrational motion of the ions is treated in Sec. IIl 



(2) 



(3) 



(4) 



Different làser cooling techniques can be used for the proper initialization of the motional state of the ions. 
They are described in Sec.^ The initial internal state where all the ions are in the state \g) can be reached by 
optical pumping to atòmic states fast decaying to the ground state \g) (Sec.0 and VI). 



The influence of the decoherence on the motional state of the ions is suppressed by làser cooling to ground 
motional states of the normal modes. The internal levels of the ions representing the qubit states |0) and |1) are 
selected such that they form slow transitions with excited states of long lifetimes. A very detailed discussion of 
the decoherence bounds of trapped atòmic ions can be found in Ref. 123] . 

Single-qubit quantum rotations can be realized on any ion and two-qubit controUed-NOT and multi-qubit 
controUed-NOT quantum gates can be applied between chosen ions duc to the poss ibility of individual addressing 
with làser beams. The implementation of quantum gates is discussed in Sec. VII, 



(5) The result of a computational process on cold trapped ions is encoded into the final state of the internal atòmic 
states. This Information can be very efficiently read out using the electron shelving method addressed in Sec. VI. 



II. ION TRAPPING 



Due to the charge of atòmic ions, we can confine them by particular arrangements of electromagnètic fields. For 
studies of ions at low energy two types of traps are used. (i) Penning trap uses a combination of static elèctric and 
magnètic fields and (ii) Paul trap confines ions by oscillating elèctric fields. Paul was awarded the Nobel Prize in 1990 
for his work on trapping partides in electromagnètic fields The operation of different ion traps is discussed in 
detall in Ref. p^ . For the purpose considered in this paper we will discuss only one trap configuration: the linear 
Paul trap (FIG^). We will follow Ref. [0 and for the mathematical treatment. 

The linear Paul trap is basically a quadrupole mass filter, which is plugged at the ends with static elèctric potentials. 
An elèctric potential 

(t,Q^Uo + VQCOs{üt) (2.1) 
oscillating with the radiofrequency is applied between two diagonally opposite rod elèctrodes. The elèctrodes are 



coupled together with capacitors so that the potential (2.1) is constant as a function of the z coordinate. The other 
two rod elèctrodes are grounded. The resulting potential at the trap axis (parallcl with the z direction) has the form 

C/o + Vo cos{üt) 



2rl 



2rl 



(2.2) 



where rp is the distance from the trap centre to the elèctrode surface. In this field the (classical) equations of motion 
for an ion of the mass m and charge q are 



qE 



or rewritten in the components 



-gV0, r = 






(2.3) 


Uq + Vo cos{nt) 


X = 


0, 


(2.4) 


Uo + Vo cos(m) 


y = 


0, 


(2.5) 




z = 


0. 


(2.6) 
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FIG. 1: Linear Paul trap in the configuration with two ring elèctrodes spa ced by 2zq. The diagonal distance between a pair of rod 
elèctrodes is 2ro. Seven ions are confined at the trap axis. The potential ( ^.l| ) is applied between two diagonally opposite rod elèctrodes. 
The other two are grounded. Ideally, equal static potential is applied on both ring elèctrodes. Used by kind pcrmission of Rainer Blatt jï^ . 



After the substitution 



b = 



2qVo 



Eq. (2.4) and ( |2.5| ) takc the form of the Mathieu equation 

+ 2òcos(2C) 



de 

dC 



nt 



= 0, 



2òcos(2C) 



= 0. 



(2.7) 

(2.8) 
(2.9) 



The Mathieu equation can be solved, in gene ral, using the Floquet solution. However, typicaUy we have a <^b^ 1, 
then the approximate stable solution of Eq. (2^) and (2^) are 



x{t) 

y{t) 



2/0 



- cos(r2í) 
^ cos(í2í) 



COS{uJyt + (fy) , 



(2.10) 
(2.11) 



where 



(2.12) 



n /ò2 n 

üJx — —\ ha, ujy — —\ a 

2V2 ^ 2V2 

and xo, yo, (px, fy are constants determined by initial conditions. We see from Eq. ( 2.10 ) and ( ^.11 ) that the motion 
of a single trapped ion in the radial direction is harmònic with the amplitude modulated with the frequency Q. 
The harmònic oscillation corresponding to the freqüències ujx and ujy is called the secular motion, whereas the small 
contribution oscillating at fi is termed the micromotion [^3| |ï^. We can eliminate the micromotion under certain 
conditions For instance, well chosen voltages on additional compensation elèctrodes (not shown in FIG.|^) nuU 
the micromotion. Then the ion behaves as if it was confined in a harmònic pseudopotential ^20 in the radial direction 
given by 

#2D = ^ + ) . (2.13) 



TypicaUy, [/q = V and hence a = 0, so the radial freqüències uJx and ujy are degenerated. Then Eq. ( 2.13 ) reduces 
to 



(2.14) 
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where the radial trapping frequency lOr is given by 

= — F = ^ 1= ■ (2-15) 

In experiments ||ï|, |ï|, |Ï5|, ^ [Ï7|, typical operating parameters are Vq ~ 300 - 800 V, n/27r ~ 16 - 18 MHz, 
ro = 1.2 mm, so we achieve the radial frequency Ur/'^tr — 1.4 — 2 MHz for Calcium ions ^°Ca^. In nature, 97% of 
Calcium consists of this isotope. To provide confinement along the z direction, static potentials Ui and U2 are apphed 
on the ring elèctrodes. Ideally, Ui = U2 = Ui2. Numerical calculations show that the potential near the trap centre 
at the trap axis is harmònic with the approximate axial trapping frequency uj^ given by 

\mulzl « í,qUi2 , (2.16) 

where zq is the distance from the trap centre to the ring elèctrode and ^ is a geomètric factor describing how much of the 
static field from the ring elèctrodes is present along the trap axis jÏ3|. Typical parameters are /27r ~ 500 — 700kHz 
for U12 — 2000 V and zq = 5 mm ||ï^, |l7|. The resulting pseudopotential for ions confined in the linear Paul trap in 
all three directions takes the form 

= ^ (x^ + J/2) + ^ (2.17) 

where the radial trapping frequency is given by Eq. ( |2.15| ) and the axial trapping frequency lo^ is defined by Eq. 



(2.16). For vàlues of experimental parameters given above, we can calculate the depth of the potential well in the 
axial direction (wz/27r ~ 700 kHz) 



? 2 
raLú,zt 



V; = ::!^~100eV (2.18) 



and in the radial direction [ujrl'í'K ~ 2 MHz) 



rrvjjzrt 



K· = ^^^~ 820eV. (2.19) 



The potential well in the radial direction is almost several times deeper than along the trap axis, i.e. there is a strong 
binding in the radial direction. Therefore we will not take into account radial oscillations of the ions in our further 
considerations. 

Finally we briefly mention how ions are loaded into the trap. We will foUow the account of practical procedures in 
Ref. [Q. Before starting the loading process, the trapping potentials are turned off for a while in order to get rid of 
any unwanted trapped residual ions. The atòmic oven producing Calcium àtoms is switched on and heats up. This 
takes about a minute. Then we turn on the electron gun ionizing neutral Calcium àtoms directly in the trapping 
volume. Cooling làsers are directed on the ion cloud containing several hundreds of ions with a diameter of about 
200 /xm. The ion cloud gradually relaxes into a steady state where the radiofrequency heating (from the elèctrodes) 
is balanced by làser cooling. The number of trapped ions is reduced by turning off the cooling. At low ion numbers, 
the ions undergo a phase transition and form a linear crystal structure. Therefore, we refer to the ion crystal or to 
the ion string or eventually to the ion chain. The loading process itself takes normally about a minute. 



III. COLLECTIVE VIBRATIONAL MOTION 



A. Equilibrium positions 

We have learnt that the ions form a linear crystal structure in the linear Paul trap after the loading process. We 
will assume a string of N trapped ions. Due to the strong binding we can neglect the radial oscillations. However, 
if a large number of ions is confined in the trap, the radial vibrations become unstable and the ions undergo a phase 
transition from a linear shape to an unstable zig-zag configuration. The relation 

ac„t=cN^ (3.1) 

determines a critical value for the ratio of the trapping freqüències a = (wz/wi-)^ for a given number of trapped ions 
iV. When a exceeds the critical value ttcrít, the ions are exposed to a zig-zag motion. The experimental vàlues of the 
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constants in Eq. (3.1) are c ~ 3.23 and (3 ~ —1.83. For experimental details and the theoretical treatment we refer 
toRef. ||. 

Further, we describe the coUective vibrational motion of the i ons. We will foUow the treatment given by James in 
Ref. The ions are exposed to the harmònic potential ( 2.17 ) duc to the trap elèctrodes and also to the repulsive 

Coulomb force from each other. Taking into account all the assumptions given above, the potential energy of N ions 
confined in the linear Paul trap is given by the expression 



N 



i=l 



N 

E 

i,j=i 



47reo \zi{t) - Zj(í)| 



(3.2) 



where Zi{t) is the position of the íth ion numbering them from left to r ight with the origin in the trap centre, m is the 
mass of the ion with the charge g, ujz is the axial trapping frequency ( 2.16 ) and eq is the permitivity of the vacuum. 
Assuming that the ions are cold enough, we can write for the position of the íth ion 



z,{t) = + A,(í), 



(3.3) 



where Zi is the equilibrium position and Aí(í) expresses small vibrations around Zi. The ions placed in the equilibrium 
positions minimize the potential energy. Hence these positions are determined by the condition 



dV_ 

dz. 



0: 



1, 



.N. 



where z = {zi, . . . , zjv) and z = (zi, 



,zn). We introduce a scaling factor 7 by the relation 
7^ 



and the dimensionless equilibrium position as Zi = Zí/7. Then one can rewrite Eq. (3.4) to the form 

i-l , JV 



E 



1 



^ {Zi ZjY j=i+i -^j) 



E 



1 



0, 



1, 



iV = 1 is a trivial case {Zi 



0). We can find the analytical solution of Eq. (3^) for two and three ions: 
N = 2, Zi = -^Ï74, Z2=^Ï74, 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Numerical calculations are necessary for > 4. For the Calcium ions ^''Ca^ and the trap frequency uüz/'ii: ~ 700 kHz, 
we may calculate the equilibrium positions as 



^ = 2, 
= 3. 



7.7 /im , 
6.1 /im . 



(3.8) 



The minimum value Azmi„ of the distance between two neighbouring ions in the trap occurs at the centre of the ion 
crystal, because the outer ions push the inner ions closer together. It has been calculated from numerical data that 
this minimum distance is given approximately by the relation pS, pO| 



Az„i„(A^) 



2.018 

^0.559 



7- 



(3.9) 



However, slightly different numerical results may be found in Ref. The relation ( |3.9[ ) happcns to be important 

when one considers individual ion addressing with a làser beam. Quantum statistics of the ion ensemble is not 
considered here because the spatial spread of the zeropoint wavefunctions of the individual ions is of the order of 
10 nm and the wavefunction overlap is then negligible |19|. 
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B. Normal modes 



The (classical) Lagrangian of the ions in the trap is given by the formula 



N N 



fc=l 



k.l=l 



dzkdzi 



(3.10) 



where we have expanded the potential energy ( |3.2[ ) in a Taylor series about the e qui hbrium positions. In the expansion 
we have omitted the constant term and the hnear term which is zero [see Eq. (3^)]. Higher order terms 0{A^) have 
been also neglected. However, they may cause a cross-couphn g bet ween difFerent vibrational modes which becomes 
a source of decoherence |2^]. The partial derivatives in Eq. ( 3.10 ) can be calculated expUcitly and we obtain the 
expression 



N 



N 



(3.11) 



,k=l 



where 



Vki = 



dzkdzi 



N 



(3.12) 



It foUows from Eq. (3.12) that Vki = Vik- The vàlues of Zj are given by Eq. ( p.7| ) for N = 2 and for = 3, whereas 
they have to be calculated numerically for iV > 4. 

The dynamics of the trapped ions is governed by the Lagrange equations 



d dL dL 



0: 



(3.13) 



dt aAfc 9Afc 

with the Lagrangian given by Eq. (3.11). We will search for a particular solution of Eq. ( |3.13 ) in the form 

Ak{t) ^ Cke-'''\ k^l,...,N, (3.14) 
where Ck are constants. Substituting Eq. ( |3.14 ) into (3.13) we get the condition for v in the form 

L·lVki-i^^SkiW^O, (3.15) 



where Ski is the Kronecker symbol and ||...|| denotes the determinant. The equation (3.15) has in general up to N 
real and nonnegative solutions Va ■ The freqüències Va are characteristic parametcrs of the system. They depend only 
on its phy sical features (not on initial conditions). A general solution of Eq. (3.13) is a superposition of particular 
solutions (3.14) and we may write 



N 



where 



By definition we will require the vectors 



a = 1, . 



to be the eigenvectors of the matrix Vki defined in Eq. (3.12), i.e. 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



k=l 
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and also to be orthogonal and properly normalized 

N 



(3.20) 



k=l 



We will number the eigenvectors in order of the increasing eigenvalues ^q. It can be shown that thc first two 
eigenvectors [a = 1,2) always have the form 



D(i) = ^(1, !,...,!), Aii = l, 
V W 

D(2) = , ^ (Zi,Z2,...,Zjv), Ai2=3. 



(3.21) 
(3.22) 



We should emphasize that Eq. (3.21) and ( 3.22| ) (they characterize two bàsic collective motional modes) are not 
dependent on the number N of the ions in the trap. Next eigenvectors (a > 3) must be, in general, calculated 
numerically. Substituting Eq. (3.21) into (3.20) we get the relation 



N 



fc=i 



(a) 



Qí = 2, , 



,N. 



We can determine analyticaUy the eigensystem for two and three ions: 



iV = 2, 



iV = 3, 



D(i) 
D(2) 
D(i) 
D(2) 
D(3) 



1 



^(1,1), m = l, 

-i=(-l,l), Ai2=3, 



^(1,1.1), 

^(-1,0.1): 

^(1,-2,1), 



M2 = 3, 
/i3 - 29/5 . 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



For larger A'^, the eigenvectors and eigenvalues must be computed numerically. The numerical vàlues for up to ten 
ions can be found in Ref. [E0|. 



Substituting Eq. ( 3.16 ) into ( 3.11 ) we get a new expression for the Lagrangian 



N 



(3.27) 



where 



(3.28) 



The Lagrangian ( 3.27 ) has split into N uncoupled terms, where Qa [Eq. ( ^.17 )] refer to the normal modes and Va 
defined in Eq. ( p.2$ ) are term ed the nor mal freqüències. Finally, the position of the ith ion in the trap can be 
rewritten in terms of Eq. (3.16D using (3.3) to the form 



Zi{t) = Zi + 



Í = 1,...,7V, 



(3.29) 



where 3fï{...} denotes the real part and Ca are constants givenby initial conditions. The collective vibrational motion 
of trapped ions determined by the eigenvector D'^' [Eq. (3.21)] refers to the normal mode called the center-of-mass 
(COM) mode 



z^^\t)^z, + n{^Cie 



í = 1,...,7V, 



(3.30) 
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and corresponds to all of the ionsoscillating back and forth as if they were a rigid body. The motion determined by 
the next eigenvector D*^^\ [Eq. ( |3.22 rcfcrs to thc breathing mode 




r7,p-'(^^V3)t > í = i,...,iV. (3.31) 



E 



N -2 
fe=l 



It corresponds to each ion oscillating with the amplitude proportional to its equilibrium distance from the trap center. 
The COM motional mode can be excited in experiments by applying an additional AC voltage on one of the ring 
elèctrodes. For exciting the breathing motional mode, a 300-times higher voltage must be applied ||ï^. Higher 
motional modes require gradient field excitation due to the nontrivial configuration of the ions in the ion string. 
However, in the limit of large ion trap dimension in comparison with the ion crystal dimension, the elèctrode elèctric 
fields are almost uniform across the ion crystal and the COM mode is very susceptible to heating due to these fields. 
Therefore, it seems to be more advantageous to use rather the breathing mode, which is much less influenced by 
uniform field s, as the quantum data bus. This will be discussed in more detail later on in the section on sympathetic 
cooling (Sec. |VQ ). On the other hand, the ions can be easily addressed with a làser beam in the COM mode, while 
higher modes require accurate bookkeeping when addressing distinct ions in the ion crystal ^ . 

C. Quantized vibrational motion 



The normal modes Qa are uncoupled in Eq. ( |3.27 ), so the corresponding canonical momentum conjugated to Qc 



is Pa — rnQa and one may write the (classical) Hamiltonian 

N N 



H-é^UPl-^'ill-lQl- (3-32) 

Oi—l a—1 

The quantum motion of thc ions can be considered by introducing the operators 



Qa ^ Qa = \l TT— («1 + àa) , (3.33) 



- Pa^^^|^^{ài~àa) (3.34) 
with the corresponding commutation relations 

[Qa,Pl3] = iflSalS , [àa, òjj] = (5q/3 ■ (3.35) 

The Hamiltonian operator associated with the external (vibrational) degrees of freedom of the trapped ions is then 
expressed as foUows {H — > Hext) 

N 

Hext = J2 («"^" + 1/2) , (3.36) 

a=l 

where àa and àj^ are the usual annihilation and crcation operators referring to the ath normal mode. We use the 
Standard notation for the number states associated with the collective vibrational motion of thc ions 

àÍàa\na) = UalUa) , (3.37) 

where \na) refers to the state of the ath normal mode and Ua denotes the number of vibrational phonons in this 
mode. The states {|«.q)} form the complete and orthonormal basis 

{ma\nfj) = ÓaP^mn- (3.38) 



We can quantize the motion of the ions by applying Eq. (3.33) to the relation (3.1(;) and expressing the displacement 
operator of the ith ion in the time-independent picture 

= EA^"^\/^(«" + ^") = E^^^^o(4 + «a), í = l,...,iV, (3.39) 

OL=\ ^ Q — 1 
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where [see Eq. {^.28\)] 



IC. 



(a) _ 



D 



zo 



2mLüz 



We can easily calculate from Eq. (3.21) that for the COM mode applies 



(1) _ 



N 



and for the breathing mode [Eq. ( 3.22| )] 



(2) 



(3.40) 



(3.41) 



(3.42) 



1=1 



Although we have not considered the radial vibrations due to the strong binding of the ions in the radial direc tion , 
a detailed treatment of the ion motion in the trap would require the extension to all three dimensions. Then Eq. (3.3) 
has to be replaced with 



qí=qí + Aqi, i = l,...,7V, 



(3.43) 



where Qí denotes the equilibrium position of the ith ion in the 3D space and Aq^ is its displacement from the 
equilibrium position. We can write 



q^ = ï^x + y^y + ZiZ, 
Aq^ = AjX + AAr+jy + A2Ar+iZ, i = \,2,...,N , 



(3.44) 



where x^, yi, zi are the equilibrium positions of the «th ion and x, y, z are unit vectors in the 3D space. The free 
Hamiltonian associated with the vibrational motion in the 3D space reads 



ilT? = E + 1/2) 



(3.45) 



and the displacement operators in Eq. ( |3.44| ) are given as foUows 



3iV 



a=l 



,3iV, 



(3.46) 



where the numerical factors /C,-"^ in general have to be determincd numcrically. 



IV. LASER-ION INTERACTIONS 



Information is encoded in infernal (atòmic) states, while it is transferred via external (motional) states of the ions. 
We can manipulate these states due to laser-ion interactions. It can be accomplished in the travelling-wave and 
standing-wave configurations. We will address in detail both approaches in what foUows. However, we should first 
comment on the selection of the two infernal atòmic levels to form the qubit. There are three possibilities ||20|| : 

• We can employ a ground and metastable fine structure excited state. This applies for ions with zero nuclear 
angular momentum [FIG.^(a)]. In this case we refer to the single beam scheme and we can drive transitions 
on optical freqüències. This configuration is used, for example, by the group in Innsbruck using Calcium ions 
40Ca+ @,|2|. 

• We can also choose two sublevels of a ground state within the hyperfine structure (ions with nonzero nuclear 
angular momentum) [FIG.^b)]. The spacing of such two sublevels is in the range of GHz. Thus, a two-bcam 
Raman scheme via a third virtual level is required in order to resolve the individual sublevels. Experiments in 
this configuration with Bcryllium ions ''Bc^ were performed in Boulder III E2l p6[. 
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• It also possible to apply a magnètic field and consider two Zeeman sublevels of the ground state [FIG.0(c)]. 
This scheme also requires Raman excitation. In this class, we can mention, for example, Magnesium ions ^"^Ig^ 
used by the group in Garching p7[ |. 

We have to mention also other active groups running experiments towards quantum lògic with trapped ions. For 
instance (in alphabetical order) IBM Almadén using -'^^^Ba"'' pl], Imperial CoUege ('*°Ca''', ^^^Hg+) p^. JPL in 
Los Angeles {^^^líg+) M, Los Alamos National Laboratory (WÍa+) H, ||l^xford University ("'"CeT) H |1, 
University of Aarhus ('^^a+) H, University of Hamburg (i38Ba+, i^^Yb^'T^ll and University of Mainz (^Ca^) 

il- 

We can use dipole and quadrupole transitions. Theoretically, the difFerence is only in the interaction constants as we 
will see later on in this section. On the other hand, in experiments quadrupole transitions have much longer lifetimes 
(one second for Calcium ions) comparing to fast decaying dipole transitions (10~*s). Experiments on an octupole 
transition in an Ytterbium ion has also been realized. The predicted theoretical lifetime in this system is of the order 
of 10^ s pTt - However, in this case one deals with very weak transitions with very stringent demands on the làser 
sources used in the experiment (although they are of major interest as potential ion trap clocks). Moreover, weak 
transitions have to be driven with a very intense làser which enhances the possibility for off-resonant excitations. 
From now on we will describe in this paper all experimental procedures for Calcium ions ^''Ca^ (FIG.0). 



co, 



le> 



(a) 




l/2F=l,mp=l 



F=2,nij=2 ' 



(b) 




1/2 inj=-l/2 



mj=+l/2 



(c) 



FIG. 2: Three possible choices for two internal atòmic states representing the qubit: (a) a ground and a metastable excited state, (b) 
sublevels of a ground state and (c) Zeeman sublevels of a ground state, where w^, lüj^i and 01^2 refer to the làser freqüències p?) . 



In the foUowing we will deal with the single beam scheme, i.e. transitions being driven by a single làser beam. 
We will not treat here the Raman scheme. The derivation of the Hamiltonian in this scheme can be found in Ref. 

. We just mention that the final Hamiltonian in the Raman scheme has the same form as the one in the single 
beam scheme, except for differences in coupling constants and for atòmic freqüències which are Stark light shifted. 
In the Raman scheme the resulting effective light ficld has the direction (frequency) determined by the difference of 
the wavevectors (freqüències) of the two participating làser beams, where each beam is represented (in a semiclassical 
approach) with a monochromatic travelling wave. Finally, the single beam scheme requires a very high làser frequency 
stability, while in the Raman scheme we only need to control the relative frequency stability between the two làser 
beams which is technically less demanding. With the Raman scheme we can also ensure the relative wavevector of 
the two beams to be parallel to the trap axis which suppresses the coupling to radial motional modes. On the other 
hand, the Raman scheme can introduce significant Stark light shifts |22[| . 




FIG. 3: Five lowest available atòmic levels of the Calcium ion ^"Ca"*' |l4[ . All transitions are accessible with sòlid state diode làsers. The 
nuclear spin of this isotope is zero, i.e. there is no hyperfine structure. For the spectroscopic notation of the levels see the text. 



In the rest of the paper we will use the Standard atòmic level notation n'^^'^^Lj, where n is the principal quantum 
number, S is the spin angular momentum, L is the orbital angular momentum and J is the total angular momentum 
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of electrons. For the fine structure case the notation is n'^^^^Lj{raj) where mj is the projection of J onto the 
quantization axis. In the case of the hyperfine structure we denote n'^^~^^Lj{F,mp) where F is the total angular 
momentum of the atom (electrons + nucleus) and mi? is the projection of F onto the quantization axis. 

Let us consider that the ion has two internal levels, denoted \g) (lower) and |e) (upper) with corresponding energies 
Eg and E^, where the transition frcquency is uq — (Ee — Eg)/h. Then the free Hamiltonian associated with the 
internal degrees of freedom is given by 

H„,t = Ee\e){e\ + Eg\g){g\ = + ^^^l.r.t , (4.1) 

where tr^ = |e)(e| — \g){g\ and íint — |e)(e| + |5)(<?|· Finally, we can write the total free Hamiltoni an fo r the jth ion 



of N ions confined in the trap communicating via one of the coUective vibrational modes [see Eq. (3.36)] 



Hoj = H,nt + Hext = ^^<7zj + hva^à , (4.2) 

where we have omitted constant terms [Ee + Eg)/2, hv /2 and dropped down the index a denoting a vibrational mode. 
The motional mode used for manipulations (especially quantum lògic operations) with the ions is callcd the quantum 
data bus because, as we will see later, it serves to transfer the Information between distinct ions within the ion crystal 
(representing a quantum register). We will consider for this purpose only the COM mode {v = oj^) or the breathing 
mode {v = ujzV^)- 

Further, we assume a powerful làser, i.e. the interaction with the ions has no influence on the làser photon 
statistics. Therefore, we will employ a semiclassical description of the làser beam. We will consider the làser beam in 
the (i) travelling-wave and (ii) standing-wave configuration. 



A. Travelling-wave configuration 



There are two different ways for addressing the ions. We can set the làser beam at a fixed position and shift the ion 
string by a very slight variation of the DC voltage on the ring elèctrode. On the other hand, we can fix the ion string 
and scan the làser across the string. In this case an acousto-optical modulator is used for làser beam deflection . 

Let us approximate the làser beam as a monochromatic travelling wave (FIG.0). We can write 



E — Eqe cos (lulí — K ■ q+ (fj 
Eoe 



(4.3) 



where Eq is the real amplitude, e is the polarization vector with |e| = 1, lvl is the làser frequency, k = kïi = (wl/c)!! 
is the wavevector with |n| = 1, c is the speed of light, q is the position vector and (j> is the phase factor. The full 
Hamiltonian for the jth ion is given by 



(4.4) 



where the interaction Hamiltonian (assuming a hydrogen-like atòmic configuration) expanded to second order (ne- 
glecting magnètic dipole interaction) has only two terms 



The elèctric dipole (DP) term is dcfined as foUows 

tDP 



(4.5) 



(4.6) 



summing over a = x,y,z. We refer to Eq. (4.6) as the dipole approximation. The elèctric quadrupole (QD) term 
reads 



QD 



a,ò 



a£:fc(í,Rj) 
dqa 



(4.7) 



where the sum is applied over a,b = x,y, z and we refer to Eq. (4.7) as the quadrupole approximation. We denote q^, 
to be the electron charge, ij is the internal position operator associated with the position of the valence electron in 
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beam 



lon 



J+1 




trap 



FIG. 4: The travelling-wave configuration corresponds to illuminating the jíth ion in the ion string with the làser beam of the frequency 
lúXj at the angle i? to the trap axis. 



the jVcí ion and Rj — (0, 0, Zj) is the external position operator corresponding to the position of the jth ion in the 
trap. 

For the present we wiU consider only the dipo le term (46), regarding to the situation when the dipole interaction 
is present and the quadrupolc contribution (4.7) is then negligible. Latcr wc wiU al so co mmcnt on the quadrupole 
interaction. If we consider consider only a single motional mode, we get from Eq. ( ^.3£ ) for the external position 
operator of the j ion 



Zj = Zj + JCjZo{à^ + à) . 

Then we can sandwich the internal position operator with the unity operator TLj — |6j)(6j| 
Eq. (4.6) to the form 



-<le 



^21 i^f,-Í'^Lt-T,jià^+à)+4'j] +H.C.| 



(4.8) 

\gj){gj\ and rewrite 
(4.9) 



where (reg)^ = {ej\rj\gj), à+j = |e ,)(gj |, à^^ = |.gj)(eJ, k = lül/c, r].j = Kj-q, 77 = k^zo, = kcosí?, 4>j ^ cj) - n^Zj 
with JCj and zq defined by Eq. (3.40). In Eq. (4.£) we consider that {ej\rj\ej) = {gj\ij\gj) = 0, because we 
assume spatial symmetry of the wavefunctions associatcd with the internal atòmic states \gj) and \ej). The schematic 
configuration is depicted in FIG.0. It is useful to transform to the interaction picture defined by the prescription 



Hm 



H = Ho + V 



n = ulv Uo 



Uo = exp 



iHot 



(4.10) 



The Hamiltonian ( ^j| ) after the transformation to the interaction picture ( 4.10| ) reads 



-o-+j exp 



-iSt 



H.c. 



(4.11) 



where 6 — ujl — ^0 and we have neglected rapidly oscillating terms at the frequency ujl + ^0 compared with low- 
frequency terms at ujl — i^o- In practice ujl — ^0, therefore to a good degree of approximation for times of interest, 
high- frequency terms average to zero [pï]| . This approximation is called the rotating wave approximation (RWA). In 
Eq. ( 4.11 ) we substitute rij = fj/y/N for the C OM mode or rjj = ??%/( V^X)/ï=i ^f) ^^r the breathing mode. The làser 
coupling constant Xj introduced in Eq. (4.11) is defined by the relation 



\DP 



qeEo 



'^{ej\{rj)a\gj)£a 



-itpj 



qeEo 



(l·eff)j 



(4.12) 



However, for a dipole forbidden transition when (ej\rj\gj) = 0, the dipole term ( |4.6| ) does not contribute (A^^ = 0) 
and the key role i s played by the weaker quadrupole interaction. In that case the làser coupling constant in the 
Hamiltonian (4.11) reads 



2hc 



^(ej|(fj)a(fi)6|.gj>ae6 



a.b 



(4.13) 
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where all parameters are defined in Eq. (4.3). 

Next, let us assume the detuning S oi the làser frequency luj^ from the atòmic frequency luq for tlie vibrational 
frequency in the form 



k = 0,±1,±2,. 



and apply the Baker-Campbell-Hausdorfï theorem |^ e^^^ = e^e^e il^^^^ Eq. ( 4.11 ). Then we can write 



(4.14) 



(4.15) 



If the làser is tuned at the frequency lül such that /c > 0, the spectral line is termed the fcth blue sideband. For k = 
the line is called the carrier and for fc < refers the fcth red sideband because the làser is red (blue) detuned from the 
atòmic frequency üjq (FIG.^. 



le>IO> 
"ig>IO> 



5 = 



le>ll> 



lg>ll> 



Ie>l2> 



Ig>l2> 



le>IO> 



lg>IO> 



(a) 



le>ll> 



lg>ll> 



(b) 



Ie>l2> 



Ig>l2> 



le>IO> 



lg>IO> 



5 = V 



le>ll> 



lg>ll> 



(C) 



Ie>l2> 



Ig>l2> 



FIG. 5: Scheme for the transition driven on the carrier (a), on the first red sideband (b) and on the first blue sideband (c), where lüIj 
denotes the làser frequency, a;o is the atòmic frequency and u is the vibrational frequency of the respective motional mode. The parameter 
S is the detuning defined as 5 = uji^ — lüq. 



When the constant Xj is sufficiently small we can assume that there are no excitations on ofï-resonant transitions 
(weak coupling regime) . Then the level structure of the ion can be considered as a series of isolated two-level systems 
pï[ . Precisely what is meant by sufhciently small is detailed in Appen dix Assuming the weak coupling regime, we 
can neglect off-resonant terms (a — /3 — fc 7^ 0) and rewrite Eq. ( 4.15 ) for fc > to the form 



n 



(+) 



^à+,(àt)i'=i^,(atà) + ^à_,^t(àtà)ai'=i 



and for fc < 



In the last relation we introduce the operator function 



{à'à 



a!(a+ |fc|)! 



(4.16) 



(4.17) 



(4.18) 



Although we allow the parameter fc to be positive or negative, we keep wri ting i ts absolute value |fc| in both cases in 
order to avoid the tricky notation of form (à^)^*^ and à^'^ in the relation ( 4.17 ) and also in what foUows next. The 
final form of the Hamiltonian is given by 



n 



(+) 



and 



n 



(-) _ 



n=0 



n.k 



\ej){gj\(g,\n+\k\){n 



n,k\ 



\gj){ej\ (g) \n){n- 



(4.19) 



|ej)(5j| ® \n){n- 



\g,){e,\<E>\n+\k\){n\ 



(4.20) 
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We have defined a new coupling constant 



n.k 



\k\ 



where 



^ ^ 7TÍ ' 



a: / n + a 
m\ \n — m 



(4.21) 



(4.22) 



is the generalized Laguerre polynomial and 
operator for time-independent Hamiltonians (4.19) and (4.20) 



n—m 



{n+ay. 
{n—m)\{a-irni)\ ' 



Úl^^ = exp - 



Finally, we may write the unitary evolution 



(4.23) 



which is given for fc > by the formula 



U. 



(+) 



= ^cosj 

n=0 

C30 

— i E^ sin 



n=0 
\k\-l 

n=0 



e,){e,\(S\n+\k\){n+\k\\ + (ffj | ® |n) (n| 



\e,){g,\ (g> \n+\k\){n\]e-^^^ + \g,){e,\ ® |„) („ + |fc| | e'^. 



(4.24) 



and for fc < 



(-) _ 



5Z '^'^^ I 

n=0 \ 

oo 

ï sin 

1*1-1 



ej){ej\(g>\n){n\ j + (^|5,)(.g,| ® |n + |fc|)(n + |fc|| 
e,){g,\® \n){n+\k\\)e-'^^^ + (\g,){e,\ ® \n+\k\){n\ ) e^^^" 



(4.25) 



n=0 



We have denoted 



1^3 2 



j\k\. For each value of k the phase factor 0, can be chosen arbitrarily for the first 



appHcation of Új^\ However, once chosen, it must be kept track of if subsequent apphcations of Új^^ are performed 
on the jth ion The real parameter |Í7J'''| is called the Rabi frequency of the transition |ej)|n) + \k\) or 

|ej)|n + |fc|) respectively. This term comes originally from the field of nuclear magnètic resonance (NMR), 

where it refers to the periòdic flipping of a nuclear spin in the magnètic field. It foUows that: 



A 47r-pulse (|fíj' '\t — Att) returns the system back to its initial state. For example 

|ej>|") ^ \ej)\n) . 
A 27r-pulse (|f2"''^|í = 2tt) changes the sign of the state. For instance 

\g^)\n + \k\)^-\g,)\n+\k\). 



A TT-pulse (|í7"''^|í tt) implies that 



\e,)\n+\k\)^\g,)\n), 



(4.26) 



(4.27) 



(4.28) 



where we set the phase factor (pj to be zero. Other cases may be easily calculated from the formulas (4.24) and 
(4.25) given above. 
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In what foUows we will assume that all motional mode s ar e in the Lamb-Dicke regime characterized by the Lamb- 
Dicke limit (Appendix^). Hence rjj introduced in Eq. (4.9) is called the Lamb-Dicke parameter. The Lamb-Dicke 
regime facilitates the ground state cooling (Sec. VB) and en ables to maintain the contrast of Rabi oscillations on a 
longer time scale [see Eq. (5^)]. Then the coupling constant ( 4.21 ) simplifics to the form 



Xj {irij) 



\k\ 



{n+\k\)\ 1 
n! |fc|! 



(4.29) 



For the purpose of coherent manipulations with internal states of cold trapped ions we will be primarily interested in 
the interaction on the carrier (fc — 0) and on the first red sideband {k — —1) which will be used for the construction 
of a wide class of quantum lògic gates. The corresponding unitary evolution operators in the Lamb-Dicke regime for 
the transition on the carrier {À) and on the first red sideband {B) may be determincd from Eq. (4.24) and (4.25) as 
foUows 



n=0 

oo 

i sin 

n=0 



2 



e.j){ej\ ® |n)(n 
ej){gj\ ® \n){n 



\9j){9j\ ® \n){n\ 

\9,){ej\®\n){n\]é^^ 



(4.30) 



and 



^co 

oo 

i sin 

n=0 



\B^\t 



ej)(ej | (g) \n){n 
ej){gj\® \n){n 



gj){gj\®\n+l){n+l\ 

\g,){e,\®\n+l){n\]e^^^ 



(4.31) 



+ l5j>(5jl® |0)(0|. 

The respective Rabi freqüències in the Lamb-Dicke limit [Eq. ( 4.29| )] for fc = and k ~ —1 are given by 

\B^\ - |A,|,7,V^. 
We could by analogy obtain evolution operators for other sideband transitions. 



(4.32) 
(4.33) 



B. Standing-wave configuration 

As an alternative approach to the laser-ion interactions we could choose a standing light field (FIG.^). One can 
place a mirror in the setup and let the làser beam reflect from it. The counter propagating waves interfere and create 
a standing-wave configuration with nodes and antinodes. However, it is experimentally very demanding to place an 
ion precisely to a node or an antinode. Let us approximate the incident làser beam as a monochromatic travelling 
wave 



Ei = E^e cos [lolí — k • q 
and the reflected beam as a counter propagating travelling wave 

Er ~ Eqs cos (wlí + K • q + i; 



(4.34) 



(4.35) 



where the reflected wave acquires an additional phase tt on the reflection at the perfect lossless mirror. Then we can 
write for the resulting standing wave 



E = E, 



E,. = 2£'oe sin (w^í + 0) sin (/t • q) 



(4.36) 



where the notation is adopted from Eq. (4.3). FoUowing Eq. ( |4.6| ) and ( [4.7] ) we can write the corresponding relations 
for the standing wave in the semiclassical representation 



E(í,Rj) = ~iEoe 



sin [Xj + Tl]{à) + a)] 



(4.37) 
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and 



= -iEoKaeb 



COS [xj +Vj{à^ + à)] , 



(4.38) 



where the new parameter Xj = K-sZj det ermin es th e position of the jth ion in thc standing wave and = (íJí/c)^^. 
The notation is adopted from Eq. (13) and ( |4.9| ). The condition Xj ~ refers to the jth ion placed in the node, 
whereas Xj = ti'/Z refers to the ion positioned in the antinode of the standing wave. 

làser 



^^^^ 




lon 



co, 



j-1 



beam 



j+1 




trap 



mirror / 

FIG. 6: The standing-wave configuration corresponds to illuminating the Jth ion in the ion crystal with the làser beam of the frequency 
üj^ at the angle i? to the trap axis. 



FoUowing the derivation for the traveUing-wave configuration, one can easily derive the Hamiltonian in the interac- 
tion picture for the standing-wave configuration. It takes thc form of the expressions ( 4.1(]| ) and (4.17), except that 
the làser couphng constant Aj and the operator function are replaced with \j and T^. They arc given in the dipole 
approximation (assuming a dipole aUowed transition) by 



-20 



^r(à^à) = e-(''?/2)sinfx, + f|fc|)^r'EH') 



a=0 



!(a + |fc|)! 



and in the quadrupole approximation (assuming a dipole forbidden transition) by 



J2{ej\i^j)a{rj)b\gj)naeb 



-i4> 



/ V OO 



2.a (à^g)" 

^' a\{a +\k\y. 



(4.39) 
(4.40) 

(4.41) 
(4.42) 



Comparing the expressions for the coupling constant in the travelling-wave configuration [Eq. ( 4.12| ) and (4.13)] 
with those ones for the standing-wave configuration [Eq. (4.39) and (4.41)], we find out that = 2|Aj'^| and 

\~Xf°\ = 2|A^^^|. The factor 2 arises from the expression of the standing wave (fl.Sq ) where we have superposed two 



travelling waves with equal amplitudes. Finally, the Hamiltonian can be written in the form given by Eq. (4.19) and 
( 4.20 ) with the coupling constant in the dipole approximation 



(ÍÍ"·'=)^^ = Af^e-KV2) 



\k\ 



{n+m 



and in the quadrupole approximation 



(4.43) 



(4.44) 



It follows from Eq. (4.43) that for the jth ion in the dipole approximation placed in the node of the standing wave 
{Xj = 0) only transitions on odd sidebands (|A:| = 2p -|- 1) are present. For the same ion in the antinode {xj = ""Z^) 
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only even sidebands {\k\ = 2p) are present, where p is an integ er or the ze ro. I n the quadrupole approximation the 
statements above are vahd in the opposite order [compare Eq. ( Í.43 ) with ( 4.44 )]. The reason for missing transitions 
in the standing-wave configuration comes from the destructive interference between the two counter propagating 
travelhng waves in the standing-wave field. 

We could easily write the coupling constant Ct"'^ in the Lamb-Dicke limit [see Eq. (4.29)]. We could also write 
the unitary evolution operator for the sta nding -wave confi guration. However, it differs from the evolution operator 
in the travelling-wave configuration [Eq. (4.24) and (4.25)] only in the coupling constant and in the phase factor, 
but it produces no fundamental problem for further applications. Thcrcfore, in what foUows we will consider the 
expressions and formulas for the travelling-wave configuration keeping in mind the way how to convert to a standing- 
wave configuration. 



V. LÀSER COOLING 



Làser cooling is the process in which the kinetic energy of àtoms is reduced through the action of one or more 
làser beams. The last decade brought ràpid progress in this research field and this effort culminated in 1997 with the 
award of the Nobel Prize in physics for làser cooling and trapping of àtoms ]^, Q . A recent review of different 
experimental techniques for làser cooling can be found in Ref. [|45| . 

One of the requirements for the practical implementation of quantum computing is the ability to prepare well defined 
initial States of the qubits ^ . In our case the qubits are represented by trapped ions with vibrational (external) and 
atòmic (infernal) degrees of freedom. Làser cooling enables the preparation of well defined initial states of motion 
and electron shelving serves for the proper initialization of the ion register. We will describe this method later on in 



Sec. VI. Làser cooling of trapped ions with the axial trapping frequency cj^ has two stages depending on the lincwidth 
F of the cooling transition |l^, Q : 

• Doppler cooling is applied when the vibrational frequency of the ions is smaller than the lincwidth of a transition 
used for cooling (oj^ < F). In othcr words, this means that the velocity of the ion duc to the trapping potential 
changes on a longer time scale than the time it takes the ion to absorb or emit a photon (strong làser driving 
is assumed). Therefore, we can assume that these processes change the momentum of the ion instantaneously. 
For < F we refer to the weak confinement regime (in the sense of weak binding of the ions to the ion trap) . 

• Sideband cooling is used for further cooling bclow the Doppler cooling limit and requires the vibrational frequency 
to be much bigger than the linewidth {uJz ^ F). Under this condition the ion develops well resolved sidebands 
and cooling to a lowest vibrational state is realized through driving a lower sideband. For ^ F we refer to the 
strong confinement regime. One can use instead a novel technique called làser cooling using electromagnetically 
induced transparency. 



A. Doppler cooling 

This stage of làser cooling is based on the Doppler effect. The technique is based on the fact that moving àtoms 
absorb photons from a counter propagating red detuned làser beam (tuned slightly below the atòmic frequency) and 
emit spontaneously in a random direction. After several such cooling cycles (absorption foUowed by spontaneous 
emission) we can write for the total momentum p of àtoms 

p^po+Y^hk^f'-^+Y^hk^f-'l (5.1) 

i 3 

where po is the initial momentum of àtoms, k^."^*'' and kj*^™^ denotes the wavevectors of the absorbed and emitted 
photons in the jih cooling cycle. We usually use a fast decaying dipole transition for the Doppler cooling, therefore 
the spontaneous emission is much faster than stimulated emission. The average total momentum of àtoms after many 
cooling cycles takes the form 

(p) =Po + ?iki(n), (5.2) 

where (n) is the average number of absorption and emission events (typically (n) ~ 10'^ — lO'') and by the definition 
ki = kj"''*^ (ki is the wavevector associated with the làser light). The spontaneous contribution averages to zero 
because it is randomly distributed over the sòlid angle 47r. If the làser is red detuned and counter propagating to the 
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motion of àtoms (po Ti ^l), then the velocity of the àtoms is significaiitly decreased [see Eq. ( p.2[ )]. For more details 
see ref. (ïï], H). 

The discussion above is vàlid for free àtoms but it also applies for trapped ions ]4^ , where the motion towards 
the làser is provided by the periòdic vibrations. The Doppler cooling limit corresponds to the final temperature 
Tdopp = fir/2fc, where F is the linewidth of the cooling transition and k is the Boltzmann constant. This temperature 
is typically of the order of mK However, the Doppler cooling limit can be also translated into the minimum 

average phonon number in the axial direction pi| 



1 



(5.3) 



where F is the natural linewidth of the cooling transition, a is determined from the angular distribution of the emitted 
radiation and S is the làser detuning from the atòmic frequency. For a dipole radiation pattern we get a = 2/5. The 
cooling is optimal for the detuning 6 = T ^ ojz- Concerning this condition we can rewrite Eq. (5.3) for a dipole 
transition to the form 



7 F 
ÏOÍZ 



(5.4) 



We have omitted the factor 1/2 corresponding to the zero point energy because it has a negligible contribution. 
The Doppler cooling limit is associated with the recoil of the àtoms at the spontaneous emission. 
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FIG. 7: Doppler cooling. The 51/2 Pl/2 the cooling transition. The cooling làser on 397 nm is red detuned by lOMHz. The ion 
spontaneously decays from the P-1/2 level. There is a 6% probability of the decay to the metastable state -D3/2, therefore the pumping 
làser on 866 nm is switched on. 



The relation (|]^) leads us to a discussion of how to choose the axial trapping frequency u)z ■ In order to be able to 
address individually each ion with a single làser beam, the minimum spacing between the ions [see Eq. ( ^3.9| )] has to 
be large enough, which requires small lOz- On the other hand, we do no t wa nt the frequency lJz to be too small in 
order to make the result of Doppler cooling as efhcient as possible [Eq. (5.4)]. Thus, the design of ion traps is also 
determined by the trade-off between these two options. 

For Calcium ions ^'^Ca.^ the S1/2 ^ P1/2 transition with the natural linewidth F = 20 MHz is used for the Doppler 
cooling. The lifetime of the P1/2 level is about 7ns and this level decays with 6% probability to the metastable -D3/2 
level (FIG.0). Therefore, optical pumping between the P1/2 and -D3/2 levels is present on 866 nm. The làser on the 
cooling transition S'1/2 ^1/2 is red detuned by F/2 ~ lOMHz. On the other hand, the pumping làser at 866nm is 
kept on the resonance in order to prevent population trapping in the superposition of the S1/2 and £'3/2 levels [Q. 
For a — 2.5, lOz/^tï ~ 700kHz and 5 — T/2 we can calculate from Eq. ( |5.3| ) the minimum average phonon number 
to be {nz)min — 3.5. This number differs from experimentally measured vàlues (which are bigger) because Eq. ( |5.3| ) 
has been derived for a two-level system while the experimental realization of Doppler cooling involves a three-level 
System. Nevertheless, {nz)min — 3.5 is still not {nz)min = required for proper operation of the quantum processor 
with cold trapped ions. Therefore a second cooling stage must be launched. 



B. Sideband cooling 



Doppler cooling represents the precooling stage in experiments with trapped ions. The final stage can be realized 
by the sideband cooling technique which may prepare the ions to the ground motional state, i.e. a well defined initial 
quantum state. Firstly, we address the bàsic idea of sideband cooling. Then we illustrate this cooling technique 
on two trapped ions in and outside the Lamb-Dicke regime. Finally, we describe how sideband cooling is realized 
experimentally. 
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In the strong confinement regime {uz ^ F) a single trapped ion exhibits in its absorption spectrum well resolved 
sidebands at lvq ± kujz {k is an integer) spaced on both sides of the carrier on the atòmic frequency lüq. Sideband 
coohng occurs when the coohng làser is tuned to a lower sideband at = ujo — kujz. In the Lamb-Dicke hmit coohng 
works efíiciently with the làser tuned on the first red sideband at = — w^- Then the ion absorbs photons of 
the energy h{ijjQ — Uz) and spontaneously cmitted photons of the average energy Hlüo — Er bring the ion back to its 
initial internal state (see Appendix ^ . In every cooling cycle (absorption + emission) the motional energy of the ion 
is damped by one vibrational quantum if hwz 3> E^. This condition iniplies that in this forní the sideband cooling 
requires the ion to be in the Lamb-Dicke limit. The whole process consists of cooling cycles in which the absorption is 
foUowed by the spontaneous emission until the ion reaches the ground motional state |n = 0) and decouples from the 
cooling làser. For the sideband cooling of the single trapped ion initially in the internal state \g) and in the motional 
state where a denotes the absorption and e stays for the spontaneous emission, we may schematically write 

\g)\n) A |e)|n-l) ^ \g)\n - l) A ... A |g)|l> A |e)|0) ^ ^0) . 

The minimum average phonon number in the axial direction that can be reached by sideband cooling is then given 
bïHI 

where the parameter a has been defined in Eq. ( ^.l[ ). It is evident that now one can achieve efficient cooling to 
the ground motional state, i. e. {nz)min — 0, assuming the strong confinement regime [lOz ^ F). The limit for the 
sideband cooling [Eq. ( ^.5D ] is constrained by the recoil of the ion and is determined by the equilibrium between 
cooling and heating processes. Heating is caused mainly by off-resonant excitations on the carrier (|5)|?t.) <-> |e)|n)) 
and on the first blue sideband (|g)|?i) |e)|rí + 1)). The sideband cooling of a single ion beyond the Lamb-Dicke limit 
also cxists and is based on the creation of a dark state in the energy level structure |Q. A single trapped Mercury 
i98jjg+ jQ-^ .^g^g gj-g^iy coolcd to the ground motional state in 1989 in Boulder while sideband cooling of a single 
Beryllium ion ^Be^ in all three dimensions was firstly reported in 1995 also by the group in Boulder |51 



1. Sideband cooling of two ions 

The key difference between one and more ions lies in the energy spectrum [ p^ pO[ . While single ions have discrete 
energy levels (three motional degrees of freedom), a chain of N oscillating ions (3iV motional degrees of freedom) 
exhibits a quasicontinuous energy spectrum due to the incommensurate freqüències of the motional modes. For 
instance, in the axial direction we have the freqüències ujz, L·lJzVS, uJzV^.S, uJzVQ-3, etc. 
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FIG. 8: Absorption spectrum /(<5) of two ions for a tlicrmal distribution P(n) in the Lamb-Dicke regime (a) and out side the Lamb -Dick e 
regime (b) as a function of detuning S (in units of the axial vibrational frequency uJz) where I{5) is defined in Eq. ( [s.fil ) and <5 in Eq. (4.14). 
Used by kind permission of Giovanna Morigi and Jürgen Eschner jj^ . 

We will discuss in detall the case of two trapped ions to illustrate the situation of sideband cooling of more than a 
single ion Q. The state of each ion (j = 1,2) will be expressed in the basis {|(7j)|n), |ej)|n)}, where n = (711,712), ni 
is the vibrational number associated with the COM mode {vi = lOz) and n2 with the breathing mode (1^2 — ^z^/ï)- 
The absorption spectrum of the jth ion will be considered in the form |55| 



E 



(n| exp (iAjKfl)|m)pP(n) 



(5.6) 
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where En = fnüzni + fiuJzVïn 2, S i s the detuning [Eq. ( [4.14 )], is defined by Eq. ( |4.9| ) and Aj is the displacement 
operator of the jth ion [Eq. ( 3.39| )]. P{n) = P(ni,n2) is a probability distribution associated with the vibrational 



motion of the ions. The Lanib-Dicke parameter distinguishes between two very difFerent regimes of sideband coohng 
of more ions: 

• In the Lamb-Dicke regime (Appendix ^ and in the strong confinement regime (w^ ^ F) only the first sidebands 
of the motional modes at wqÍWz and ujqÍuJzV^ appear around the significant carrier peak at in the absorption 
spectrum [FIG.^(a)]. The higher sidebands are suppressed due to their strength being proportional to higher 
powers in the Lamb-Dicke parameter than r/ denoted as 0{r]^). Tuning the làser on the first red sideband of the 
COM mode (6 = —ujz) we can reach its ground state \ni — 0) at the same coohng rate as for a single ion [ p4| . 
However, in the case of two ions the breathing mode is decoupled from the COM mode and its coohng is almost 
frozen. Simultaneous cooling of more modes requires the modes to be coupled to the cooling làser and then 
the requirement of the strong coupling regime {ujz ^ T) has to be reconsidered or one has to use alternative 
techniques. 

• Outside the Lamb-Dicke regime higher sidebands with the strength proportional to 0{rf') also contribute and 
the absorption spectrum exhibits the structure with many overlapping sidebands at cjq i kuz ± lujzVS where k, l 
are integers [FIG.^(b)]. In this situation the làser tuned on a lower sideband (it does not have to be strictly the 
first red sideband of the COM mode) excites simultaneously all sideband transitions around this lower sideband 
in the interval of the linewidth V. Then the COM and the breathing mode are coupled and cooled at once. 
However, the cooling process is much slower in comparison to cooling of a single ion beyond the Lamb-Dicke limit. 
It is partly caused by (i) the increasement of the number of the motional modes but also by (ii) the appearance 
of dark states . The dark states are almost decoupled from a resonantly excited state because their motional 
wave function after the absorption overlaps with the motional wavefunction of the excited state only a little. 
Thus, the ions may be trapped in these dark states and it slows the cooling process down. This problem can 
be solved by escaping from the strong confinement regime (ojz 3> F), i.e. by increasing the linewidth F. It will 
cause that a single level will be coupled to more levels (more sidebands are in the resonance) and the ions will 
be cooled more efficiently due to more cooling channels. As a result the dark states will disappear because more 
channels provide more ways for the ion to escape from dark (population trapping) states. Moreover, the rate of 
the cooling cycles (absorption -f emission) is proportional to the linewidth F. Summarizing both eífects we can 
conclude that the total cooling time beyond the Lamb-Dicke limit can be shortened significantly for F ~ a;^ . 



2. Experimental sideband cooling 



Two ions were cooled for the first time to the ground motional state in 1998 in Boulder. It was achieved on 
Beryllium ions ^Be^ illuminating both ions at once |Q. However, it is sufficient to illuminate only one ion from 
the entire ion string because other ions are cooled sympathetically due to the strong Coulomb coupling. Although 
we need only one motional mode (COM or breathing axial mode) as the quantum data bus, which has to be in the 
ground motional state, we require also other modes to be cooled close to the ground state. Uncooled motional modes 
with thermal phonon distributions significantly affect the Rabi frequency in the data mode and spoil the fidelity of 
the coherent state manipulation [sec Eq. (5.S)]. 




FIG. 9: Sideband cooling with cooling and pumping transitions. The blue transitions correspond to the cooling cycle. The làser at 
729 nm is tuned on the first red sideband. The làser at 854 nm couples the metastable -D5/2 level with the fast decaying P3/2 level. The 
cooling cycle is closed by the spontaneous emission on the fast S1/2 ^ Pi/2 dipole transition. The P3/2 level may decay to the state 
The ion is recycled to the state P1/2 by the làser at 866 nm foUowed by the spontaneous emission back to the state S'i/2- Decay to 
the sublevels 5i/2("^=+l/2) is counteracted by driving the ^ Pi/2 transition with the a~ polarized làser at 397 nm. 
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The group in Innsbruck has reahzed difFerent approaches in sideband coohng of two Calcium ions [g^ y_7|. If they 
cool only one motional mode, while the other modes are left in the thermal states, they achieve the ground state 
population greater than 95% ((n) ~ 0.05) in the respective mode. However, they can cool sequentially all motional 
modes close to the ground state. For this purpose they use a small modification in the sideband cooling scheme. The 
làser frequency and làser power has to be set sequentially for the respective first red sideband of the given motional 
mode. After sequential cooling all modes the corresponding average phonon numbers are from (n) ~ 0.05 to (n) ~ 2.3 
because the recoil energy from the spontaneous emission in the cooling process of one motional mode reheat other 
modes. 

For Calcium ions '^''Ca^ the quadrupole transition between the two Zeeman sublevels Si/2{m,j = —1/2) and 
Dç^/2{'mj = —5/2) is used for the sideband cooling (FIG.^. The làser on 729 nm is tuned on the first red sideband 
of the respective motional mode and a weak magnètic field is applied for Zeeman splitting of energy levels. The 
lifctime of the metastable D5 /2 level is about a second, therefore there is the pumping on 854 nm to the fast decaying 
^3/2(™J — —3/2) level in order to decrease the duration of one cooling cycle, i.e. to increase the cooling rate. The 
P3/2(mj = —3/2) level decays spontaneously to the initial state Si/2{mj = —1/2) and closes the cooling cycle. 
However, the P3/2 level may decay with a small probability to the -D3/2 level, therefore the pumping làser on 866 nm 
recycles the population to the P1/2 state which decays to the S'1/2 state. The S1/2 <-> P1/2 transition is driven with 
the cr^ polarized làser on 397 nm to counteract the population of the Si/2{mj = +1/2) level [ïi) . 



C. Sympathetic cooling 



In the previous section we have mentioned that it is sufíicient to illuminate with cooling làsers only one ion from 
the ion string because the other ions are cooled sympathetically due to the Coulomb interaction between them. Hence 
the term sympathetic coolin g. H owever, instead of identical ions one can consider different atòmic species (eventually 
isotopes) in the ion crystal [p7| . Then the addressing of cooling ions avoids the disturbance of infernal states of lògic 
ions which store the Information l5q| . 
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FIG. 10: Normalized axial freqüències ( 3.28 ] as a function of Ç for (a) three, (b) five, (c) seven and (d) nine ions. Used by kind permission 
of David Wineland and David Kielpinski |59| . 



Electric fields from the trap elèctrodes are one of the sources of the motional decoherence of the ion crystal due 
to heating of coUective vibrational motional modes (normal modes). If one assumes the dimension of the ion trap to 
be much larger than the dimension of the ion crystal, then we can expect the elèctrode elèctric fields to be nearly 
uniform across the ion crystal. Such uniform fields influence and heat only coUective motional modes involving the 
centre-of-mass (COM) motion of the ion crystal. Uniform elèctric fields can directly heat up the normal mode used 
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for quantum lògic as the quantum data bus. We can overcome this constrain by selecting a specific normal mode for 
quantum lògic which is decoupled from heating. However, not all motional modes are prevented from heating. In 
what foUows we will discuss this point foUowing Ref. pof . 

Let us consider the ion crystal with an odd number N of ions which consists of — 1 ions of mass m and of a central 
ion of the mass M dcfining the ratio Ç — M/ra. Now we can follow the lines in Sec. II] and find the normal modes and 
freqüències of the ion string with unequal ions. We find out that (i) there are [N — l)/2 axial normal modes for which 
the central ion does not move and corresponding eigenvectors D^") and eigenfrequencies Va = í^z^JV^ do not depend 
on the parameter C^. Moreover, these modes do not have a component associated with the axial COM motion. (ii) 
There are also other (TV + l)/2 modes having a component of the COM motion and coupling to any uniform elèctric 
field which causes their heating. 

For very small or very large vàlues of C the motional modes become degenerate and pair up (FIG.|ÏO|). From this 
point of view the value ^ ~ 1 seems to be suitable. It has also been calculated that those modes having the central 
ion at rest (neglecting gradient elèctric fields) do not heat at all (we refer to these as cold modes) , while all other 
motional modes heat to some extent depending on the value of Q. Their heating rate (average number of phonons 
gained per second) drops rapidly for C ^ 1 (FIG.|ÏÏ]). We refer to these modes as hot modes. 

Thus, it seems that the optimal choice is C ~ 1. That means the central (cooling) ion should be chosen such that 
it is identical to — 1 other (lògic) ions or is an isotope of lògic ions. For C ~ 1 we can choose the lowest cold mode 
(the second lowest motional mode called breathing mode) to be used for quantum lògic as the quantum data bus 
because in the case of C — 1 only the lowest motional mode (corresponding to the COM mode for equal ions) will 
heat significantly and can be cooled via the central cooling ion. If the value of C, differs very much from 1, we have to 
cool all [N + l)/2 hot modes via the central ion. 
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FIG. 11: Normalized heating rates for hot axial modes as a function of Ç for (a) three, (b) five, (c) seven and (d) nine ions. Used by kind 
permission of David Wineland and David Kielpinski [59[ . 

The group in Garching runs experiments where the ion string consists of Indium (^^^In^) and Magnesium (^^Mg^) 
ions. The numerical analysis for the ion crystal containing these two atòmic species ordered in different configuration 
can be found in Ref. [S^ . The mass ratio is C = 4.6 which is not within the optimal range discussed above (^ ~ 1). 
On the other hand, it can quite advantageous because the heavy ion fulfiUs the Lamb-Dicke limit (Appendix ^) 
easier than the light ion. Distinct atòmic species can also have a very different atòmic spectrum which may be found 
convenient when làser addressing closely spaced ions. However, for the heavy central cooling ion we pay the price in 
the form of hea ting rates of higher motional modes. By all means, Indium ions can be efíiciently cooled to the ground 
motional state []60| and Magnesium ions can serve for quantum lògic operations and storing the Information. 

Finally, we have to mention that the demonstration of sympathetic cooling using two different atòmic species is 
very demanding on current experimental technology due to problems of loading the ion trap with distinct àtoms in 
a desired configuration. 
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D. Làser cooling using electromagnetically induced transparency 



Quantum computing with cold trapped ions requires one of the motional modes (the one used as the quantum data 
bus) to be cooled to the motional ground state and other modes to be inside the Lamb-Dicke regime. For this purpose 
one could eventuaUy use Dopple r co ohng assuming the axial trapping frequency comparable with the line widt h 
of the coohng transition T [Eq. (^.4D]. However, it would cause very close spacing of the ions in the trap [Eq. (^)] 
with difhcultics at individual addressing with the làser beam and optical resolving. On the other hand, we can use 
sequential sideband cooling of the motional modes described in Sec. VB. However, cooling of one motional mode 
causes heating of the other modes. Moreover, sideband cooling requires a very narrow bandwidth to excite the first 
red sideband of the respective motional mode only. Otherwise, off-resonant transitions (especially carriers) are also 
driven what causes heating as well [iïs!. 
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FIG. 12: (a) Levels and transitions of tiie cooling technique using electromagnetically induced transparency (EIT). The inset sliows the 
absorption on \g) <-> |e) while strongly driving the transition \r) ^ |e). (b) Absorption on the carrier (|n) — > |n)) and on the first sidebands 
(|n) — > |n it 1)). Used by kind permission of Giovanna Morigi and Jürgen Eschncr teli. 



A novel cooling technique was developed in 2000 with a lower cooling limit than Doppler cooling and with a wider 
cooling bandwidth than sideband cooling. It was named làser cooling using electromagnetically induced transparency 
(EIT) |§, ~ 
resonance 



excited state |e) [FIG. 
Ar- The transition \g) 



32[. It is based on a quantum interference effect called EIT or coherent population trapping or also dark 
631 . It employes a three-level system with a ground state \g), a stable or metastable state |r) and an 
~2|(a)]. The transition |r) <-> |e) is driven with a làser beam of the intensity Ir blue detuned by 
í-> |e) is coupled by a weak làser with the intensity Ig (where Ir/Ig — 100) also blue detuned 
ilr\'^ oc Ir introduces a significant Stark light shift Aw where |Q 



by Ag. The intense làser with 



\Ar\)/2. 



(5.7) 



Thus, the làser on the transition 



|e) designs the absorption spectrum seen by the weak làser on |r) |e) via 



the level |e). Then there is a broad resonance at An ~ 0, a dark resonance (EIT) at Ag = A^ and a bright narrow 



resonance at A,. 



Auj ísee the inset in FIG. 



of freedom, (ii) setting the detunings such that A, 



t2|(a)]. 

Ar 



Therefore, (i) taking into account also the motional degrees 
and (iii) setting the Stark light shift equal to the vibrational 



frequency (Aw ~ z/) we obtain the absorption spectrum depicted in FIG.|ï^(b). We see that the absorption on the first 
red sideband (cooling transition) is enhanced while the absorption on the carrier (heating transition) is eliminated. 

The bright resonance width can be wide enough to cover several motional modes which can be consecutively cooled 
at once. It was experimentally demonstrated in Innsbruck on two motional modes separated in the frequency by 
1.73 MHz. The modes were cooled to their ground motional states with 74% {{n) ~ 0.35) and 58% ((n) ~ 0.72) 
occupation | |6^ . A great irnprovement of these results should possible for a rebuilt apparatus allowing optimal access 
of làser beams to the ions . 

FoUowing the advantages of làser cooling using EIT it has been estimated that all 3A^ motional modes can be 
cooled to a mean phonon number (n) < 1 for w^/27r = 700 kHz and = 10 j6^. It is very important to cool 
3A^ — 1 spectator motional modes to the Lamb-Dicke regime. Otherwise, thermally excited spectator modes cause the 
fractional fluctuations (blurring) in the Rabi frequency of the mode used as the quantum data bus for quantum lògic 
operations p^. These fluctuations in the Rabi frequency of the ath mode can be estimated as p6| 



Aü 



/E 4 ("z^) ( + 1) 



(5.8) 
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where rj/j is the Lamb-Dicke parameter of the /3th motional mode and (np) is the respective average phonon number. 
The ratio (Í7/AÍ7) determines the maximal number of Rabi cycles |ï^. A detailed description of experimental làser 
cooling using EIT on Calcium ions can be found in Ref. . 

VI. ELECTRON SHELVING 

Electron shelving is the experimental method for the discrimination between two electrònic levels with an efficiency 
approaching 100%. It was firstly demonstrated in 1986 |6^. Let us assume a three-level atom consisting of a ground 
level \g), a metastable excited state |e) and an auxiliary excited fast decaying state \r) (FIG.|Ï^). The |g) <-> |e) 




FIG. 13: Electron shelving. The \g) <-> \e) transition is coupled with a làser pulse forming the superposition a\g) + 0\e). The \g) <-> |r) 
transition is driven with a strong làser. The fluorescence detection signal is coUected on the \g) ^ |r) transition if the ion coUapses to the 
state \g). If the ion is shelved in the dark state |e), no fluorescence is observed bccause the |e) state is a metastable state. 

transition is coupled by a weak làser forming a superposition a\g) + /3|e), while the \g) ^ \r) transition is driven with 
a strong làser. If the atom coUapses to the \g) state during the measurement, a strong fluorescence signal is collected 
on the fast transition \g) <-> \r), i.e. the atom is excited from \g) to |r) and spontaneously decays back to the \g) 
state what is observed as the fluorescence. However, if the atom stays shelved in the metastable excited state |e), no 
fluorescence can be observed on the driven \g) <-!■ |r) transition. Hence the name electron shelving. Even though the 
detection efficiency is low, we can keep exciting the measuring transition \g) <-> |r) and detect some spontaneously 
emitted photons. Thus, we are able to discriminate the \g) and |e) states with almost 100% efficiency. We can also 
obtain the occupation probability |ap for the \g) state and |/3p for the |e) state averaging over many repetitions of 
the same experiment [nsl. 

In the casc of Calcium ions, the ground state \g = 5*1/2) and and the metastable excited state |e = D·^^j-i) form the 
qubit Jl^ . The auxiliary state corresponds to the P1/2 level. The 5*1/2 *^ 1^5/2 transition is illuminated with a weak 
làser pulse on 729 nm and the 5i/2 -P1/2 transition is driven with a strong làser at 397 nm. However, the ion can 
decay from the P1/2 level with a small probability to the level. Therefore, there is pumping on the P1/2 <-» £^3/2 
transition at 866 nm (see FIG.0). 

To clarify the efficiency of the electron shelving method we report briefly some results measured in Innsbruck [ p^ . 
When the ion is found in the 5i/2 state it scatters about 2000 photons in 100 ms to the detector. However, for the ion 
in the dark I?5/2 state the number of events drops to only about 150 photons in 100 ms. These 150 photons appear 
duc to dark counts of the photomultiplier and some scattered light from the làser at 397 nm. The ion string in the 
linear ion trap may represent a quantum register, where the interna! state of each ion, i.e. the state of the qubit, can 
be detected using a CCD camera. Then the ion in the \g = S1/2) state appears as a bright spot or a dark spot if the 
lon IS found in the |e state ||ï|, |ï|, |ï| . 

VII. QUANTUM GATES 

One of the requirements for the physical implementation of quantum computing in a ccrtain quantum system is a 
set of quantum gates that can be realized in the quantum system under consideration. It has been shown that any 
unitary operation can be composed of single-qubit rotations and two-qubit controUed-NOT gates . In what foUows 
we will describe how these and some more complex quantum gates can be implemented on cold trapped ions. We will 
use the notation \g) and |e) of the logical states for the qubit rather than |0) and |1) due to the representation of the 
qubit by the interna! states of the ion. 
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A. Single-qubit rotations 



A general single-qubit gate corresponds to a unitary evolution operator that acts on a single qubit and is represented 
in the basis {jg), |e)} by the matrix 



, Wgg Wge 

' Weg Wee 



(7.1) 



A special case of the single-qubit gates is a single-qubit rotation (FIG.jïJ). Its parameterization depends on the choice 
of the coordinates on the Bloch sphere. We will dcfine it in the matrix form in the basis {|.g), |e)} as foUows 



T^99 




-i 




TZee ) 





cos(6'/2) e''^sin(6'/2) \ 
-e-''^sin(6l/2) cos(6'/2) ) ' 



(7.2) 



where 9 refers to the rotation and (p to the relative phase shift of the states \g) and |e) in the corresponding Hilbert 
space. 



R 



FIG. 14: Schematical representation of a single-qubit rotation. R is defined by Eq. ( |7.2| ) in the basis {I3), |e)}. 

The single-qubit rotation can be perf orme d on a selected ion from the ion string in the Lamb-Dicke regime by 
applying the unitary evolution operator (4.3C). We may rewrite this operator to the form 



A 



= ^cos(£7r/2) 

00 

+ ^sin(^7r/2) 



',){e,\®\n){n\ j + y\gj){gj\® \n){n\ 

j){9A ^ \n){n\)e-^'f'^ + (\g,){e,\ ® \n){n\ je'^^' 



(7.3) 



where we have applied the arbitrary choice of the phase factor — > -f 7r/2) with respect to the remark below 
Eq. ( 4.251 ). The operator (7.3) corresponds to the jth ion illuminated with the làser beam on the carrier (cj^ — loq) 
with the làser pulse duration t = £7r/|Aj|, where the làser coupling constant Aj depends on the type of (i) the driven 
transition and (ii) the làser configuration (see Sec. IV). We will refer to the operation expressed by Eq. (7.3) as the 
^TT-pulse on the carrier. 



B. Two-qubit controlled-NOT gates 

A two-qubit controlled-NOT (CNOT or XOR) gate acts on two qubits denoted as a control and a target qubit 
(FIG.|Ï^). If the control qubit (mi) is in the state |e), then the state of the target qubit (7712) is flipped. Otherwise, 
the gate acts trivially, i.e. as the unity operator 1. We may characterize this gate with the help of the foUowing truth 
table 



\9r,n)\9,n2) ' \97m)\9n 

\9mi)\Sm2) * |5'ni)|Cn 

|^mi)|í/m2) ^ 1^7711)1^77 

|677ii)|Cm2) * |677ii)|5" 



(7.4) 



The implementation of the two-qubit CNOT gate on two selected ions in the ion string requires the introduction 
of a third auxiliary infernal level \r). In the original proposal the selective excitation of two sublevels of the |e) 
level is used instead. The selection depends on the làser polarization, where the \e,p = 0) and \e,p = 1) sublevels are 
considered. There have also appeared proposals how to avoid the establishment of the auxiliary infernal level to the 
scheme [^8|, ^ ^ . Two of them will be discussed later on in his section. 

Now we are ready to write two unitary evolution ope rator s corresponding to làser pulses driven on the first red 



sideband in the Lamb-Dicke regime on the jth ion [Eq. (4.31)] between the infernal levels \g) ^ \e) with the atòmic 
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m. 



m-, 



FIG. 15: Schematical representati on o f a two-qubit controUed-NOT (CNOT) quantum gate. The mi (r?i2) qubit is control (target). The 
gate is defined by the truth table ( [i'.4| ). 

frequency (jJq^ — [E^ ~ Eg)/h and for \g) ^ \r) with cjg^ — {E^ — Eg)/h. They are given by 



= cos(€^/2) 



\e,){ej\ (g> 



+ l5,>fel®|l)(l| 



- ïSÍn(£7r/2) 



e,){gj\®\0){l\ e 



|g,)(e,|®|l)(0| e'*. 



\g,){g,\®\0){0\ + 



and 



- cos(Í7r/2) 



+ |5,>(5,I®|1)(1| 



isin(£7r/2) 



kj)fel®|0)(l| e 



|5,)(r,|®|l)(0| e'*- 



+ \gj){g,\®\0){0\+0, 



(7.5) 



(7.6) 



wherewe have applied again the arbitrary choic e of the phase factor {(f)j = (jjj —Tr/2 (j)j). The symbol O in Eq. ( |7.5| ) 
and ( |7.6| ) correspond to the terms in Eq. (4.31) associated with the dynamics on higher vibrational levels for n>2. 
We do not have to consider them because the ions are assumed to be cooled to the ground motional state \n = 0). 
We use the Hilbert space spanned only by the motional states |n = 0) and \n = 1) forming an auxihary qubit used as 
the quantum dat a bu s. 



The operators (7^) and (7^) correspond to fcTr-pulses on the first red sideband (w^ 



V and ujl 



for n = with the làser pulse duration t = £TT/\Xj\T]j. Finally, the two-qubit CNOT gate on two ions corresponds to 
the evolution operator sequence (acting from right to Icft) ||] 



■^m2 \" ) "mi '-'m2 *-^mi •^m2 > ' 



(7.7) 



where À^^i^) ^-iid ■Miti'^) are given by Eq. ( |7.3| ) and stand for the 7r/2-pulses on the carrier {lül — í^q^) on the TO2th 
ion with the phase = and — tt, respectively. The operator B}^^ is defined by Eq. (7.5) and represents the 

Uz) on the mith ion with the phase factor — 0. The operator E^^J 



íl/2/ 



TT-pulse on the first red sideband (wl 

defined in Eq. (^) stands for the 27r-pulse on the first red sideband {ujl 
The middle sequence B}^^^ the evolution operator ( \l-T\ j can be schematically represented as foUows 



on the TO2th ion with 



0. 



l5mi)|5m2)|0) 

|5mi)|e„i2)|0> 

|emi)|S'm2)|0> 

|emi)|e,„2)|0) 



I5mi)l5m2)|0) 
|gmi)|em2)|0> 
-ï|3mi)|gm2)|l> 
-í|.9mi)|em2)|l) 



1712 



l5mi)|5m2)|0) 
|ffmi)|e™2)|0) 
ï|5mi)|5m2)|l) 

-ï|.9mi)|e,„2)|l) 



I.9mi)l5m2)|0) , 
Iffmi)|em2)|0) , 
|emi)|5m2)|0) , 

|emi)|e„,,)|0) . 



(7.8) 



Finally, the evolution operator (7.7) refers to the transformation 



Iffmi>l5m2)|0) 
I5mi)|em2)|0) 
|emi)|£'m2)|0) 

|emi)|e™J|0) 



I5mi)l3m2)|0) 

Iff™i)|em2)|0) 

|emi)|em2>|0) 
|emi)|gm2)|0) 



(7.9) 



on two selected ions labelled as mi and m2 in the string of N ions. 
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It is evident from the discussion above that this realization of the CNOT lògic gate on the ion system requires 
the ions to be cooled to the ground motional state |n = 0) in order t o m aintain the fidehty of the computational 
process. Otherwise, as the ions heat up, higher terms O in Eq. (^) and (7^) also contribute and introduce significant 
imperfections into the implementation of the quantum gate. The two-qubit CNOT was firstly demonstrated in Boulder 
in 1995 ||5Ï|. A singlo Bcryllium ion was uscd, where the control qubit was stored into two lowest vibrational states 
\n = 0) and \n = 1) and the target qubit was represented by two hyperfine levels \g = Si/2{F ~ 2, mi? = 2)) and 
|e = 5i/2(F=l,mi. = l)). 



C. Alternative implementation of two-qubit controlled-NOT gates 

1. Simplified quantum lògic 

Monroe et al. have proposed the realization of the two-qubit quantum lògic gate based on the precise setting of the 
Lamb-Dicke parameter [|9|. The control qubit is assumed to be encoded into two lowest vibrational states \n — 0) 
and \n = 1) of the considered collective vibrational mode, while the target qubit is represented by two infernal levels 
\gj) and \ej) of the jth ion from the string of N ions in the linear Paul trap. The CNOT gate under consideration is 
then described by the truth table 



|0)|ff,) 
|0)|e,) 



|0)|ff,), 
|0)|e,), 
|l>|e,), 
ll)!^.)- 



(7.10) 



Further, we adopt the main idea of the original proposal |69(|. Driving the jth ion with the làser on the carrier 



[u!]^ = ujq) is described by the evolution operator (4.24). The coupling constant fíj'*^ is introduced by the expression 



(4.21). Then for n = and n = 1 with fc = we get 



Let us set the Lamb-Dicke parameter such that 



(7.11) 
(7.12) 



(7.13) 



where p is an integer. The realization of the transformation (7.10) requires driving the carrier transition on the jth 
ion with the duration t such that 



We can calculate using Eq. (7.13) that 



(7.14) 



(2p- l)7r. 



(7.15) 



Substituting Eq. (7.14) and (7.1Í ) into Eq. (4.24) we find out that the infernal state of the jth ion is flipped only if 



the collective vibrational state is \n = 1). We can write 



|0)|ff,) |0)|g,) 



|1>|5.) 



IDIe 



+ÍC, 



(7.16) 



'^1l>l5,> 



This transformation corresponds to the CNOT gate (7.10) apart from the phase factor which can be eliminated 
by the appropriate phase settings of subsequent operations. 

The CNOT gate between distinct ions representing two lògic qubits can be implemented (using the proposal being 
discussed) by two additional làser pulses on the first red sideband {ul — ujq — i^)- We have on mind the CNOT gate 
given by the truth table ([7.9|). Firstly, we apply a 7r-pulse on the first red sideband on the mith ion (corresponding 
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to the evolution operator ( 4.25| ) with íl^'^'^í = tt) mapping the internal state of this ion onto the coUective vibrational 
state. We can write 

l5mi>|0) - 

|emj|0) - 

l5mj|l) - 



l.9"M>|0), 



(7.17) 



-íe-**"M|e,„J|0) 



Secondly, we apply a làser pulse on the carrier on the m2th ion representing the reduced CNOT gate (7.10) and finally, 
we map back the coUective vibrational state onto the internal state of the mith ion by reapplying a 7r-pul se o n the first 
red sideband on this ion. This sequence of three làser pulses corresponds to the complete CNOT gate ( [7.9[ ) between 
two distinct ions with the appropriate choice of the phase factors. 



Comparing this scheme to the original proposal of Cirac and ZoUer (7.7), we need fewer làser pulses to realize 
a two-qubit CNOT gates on trapped ions and there is no need for a third internal auxiliary level. However, more 
important is the overall time needed to complete the gate and the sensitivity to imprecisions. The main limitation 
of this Monroe scheme is that it is slow compared with other methods at the same level of infidelity (caused by 
off-resonant transitions) and it is rather sensitive to imprecision in the làser intensity (|ííp oc /) 



2. Fast quantum gates 

Jonathan et al. have proposed another alternative realization of two-qubit quantum gates on cold trapped ions 
scalable on N ions [0. It is based on (i) using carrier transitions and (ii) taking into account Stark light shifts 
of atòmic levels. FoUowing Ref. |^ the bàsic idea of fast quantum gates is that the resonant driving of a carrier 
transition (|5)|íi) ^ |e)|".)) with an intense làser causes the splitting of dressed states |±) = l/^/2{\g) ± |e)) in the 
interaction picture by amount 2híl, where the coupling constant Í7 is proportional to the làser intensity / [FIG.p^(a)]. 
When we set the làser intensity such that the splitting of the dressed states |— ) and 1+) is equal to one motional 



le>. 



2Í2 



l+> 



l-> 



l+>IO> 



l->IO> 



(a) 



(b) 



l+>ll> 



-1 l->ll> 



FIG. 16: Fast quantum gates on cold trapped ions are based: (a) on the splitting of the dressed states |it) = l/\/2{\g) it |e)) by amount 
proportional to the làser intensity [Q = when the làser is tuned on the carrier transition {w]^ = a;o) and (b) on setting the splitting 

such that it is proportional to one motional quantum {u = 2Q). 



quantum Hi'. Then Rabi oscillations appear between the state |+)|0) and |— )|1) with |0) and |1) referring to the 
lowest coUective vibrational states of the ions [FIG. 0(b)]. Using this swapping between the |-l-)|0) and |— )|1) states 
one can construct a CNOT gate between two distinct ions foUowing the truth table ( 7.9). Quantum gates using this 
idea are faster than Standard quantum gates on trapped ions (discussed in Sec. VII B ) approximately by the factor of 
l/ry assuming the Lamb-Dicke regime. The speed of quantum gates will be discussed in Sec. IX. 



D. Multi-qubit controlled-NOT gates 

A multi-qubit controUcd-NOT gate is defined by analogy to the two-qubit CNOT gate. The only difference is the 
number of control qubits (FIG.[Ï7|). The multi-qubit (controlled)''-NOT gate acts on q-|-l qubits with q control qubits 
(mi, . . . , ruq) and the m^+ilh qubit is target. If all control qubits are in the state |e), then the state of the target qubits 
is flipped. Otherwise, the gate acts as the unity operator 1. The truth table of the multi-qubit (controlled)'-NOT 
gate acting on mi, . . . , m^+i qubits is 



|*no)|gm,+ i) 
|*no)|em,+ i) 
|*í,eí)|ffm,+ i) 



I í't/es) |em<j_|_J 



\^no)\grn^+i) , 
|*«o)|e™,+ i) , 
|í'aes)|em,+ i) , 

|*aes)|5m,+ i) ■ 



l*a-> = n ®|em,), 



(7.18) 
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The multi-qubit (controlled)'^-NOT gate acting on g + 1 ions (toi, . . . , níq ions represent the control qubits, while the 
77iq_|_ith ion stands for the target qubit) can be reahzed by applying the evolution operator (acting from right to Icft) 















_J=2 







(7.19) 



where the B oper ators are taken for the value = 0. However, this choice of the phase factor has no fundamental 
importance. Eq. (7.19) apphes for three and more ions and the scheme requires again the auxihary qubit encoded 
into two lowest levels |n = 0) and \n ~ 1) of the coUcctive vibrational mode used as the quantum data bus. 



m, 



m„ 



m 



FIG. 17: Schematical representation of a multi-bit (controllcd)'-NOT gate acting on q + 1 qubits with q control qubits and mqj^ith qubit 
to be target. The gate is defined by the truth table (7.1Í). 



Now we verify whether the evolution operator ( 7.19 ) corresponds to the truth table of the multi-qubit CNOT gate 
given by Eq. ( 7.18 ). At first we consider only the B operators and then we comment on the action of the A operators. 
It holds for q + 1 ions involved in the multi-qubit CNOT gate that: 

• If the mith ion is in the ground state \gmi), then the action of the B operators in Eq. ( 7.19| ) corresponds to the 
unity operator. 

• If the rriith ion is excited with all other ions in the ground state \erm)\g)'^ , we get 



(7.20) 



Thus, the transformation is performed on the mith ion and then on the first next ion in the ground state. The 
state of all other ions in the ground state is not transformed. If more ions (besides the mith one) are excited 
(except if they all are excited) , their state does not change because the B^^^^ operator acts only in the Hilbert 
space spanned by {|g„i.), jr^^)} [sec Eq. (Q]. 

• If all the ions are excited, i.e. |e)'?+^, it foUows that 



)«+i|0)^-*|5„H)|e)11)- 



)«+i|0) 



(7.21) 



Finally, the A ope rato rs complete the operation (7.19) such that it corresponds to the transformation (7.18) by analogy 
to Eq. (ItJ) and (ItJ). 



E. Multi-qubit controlled-iï gates 

A multi-qubit (controlled)'^-i? gate acts again on g + 1 qubits. However, it performs a single-qubit operation ([7^ 
on the mg4_ith (target) qubit if all rni, . . . , ruq control qubits are in the state |e). Otherwise, it acts trivially (FIG. 
Speaking precisely, if all control qubits are in the state |e), then the rotation R = K[ a FÚ2 cr R2 Ri is applied (from 
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right to left) on the target qubit. In the basis of the target qubit {\g)m +i, |e)m we introduce the matrices 



R = 



cos 9 e*^*^ sin í 
-e-'^'/'siííd cos 61 



1 

1 



Ri = 



e"^ 
-e"*-^ 



oÍ0 



(7.22) 



Ra = 



cos(6l/2) sin(6'/2) \ 
-sin(6l/2) cos(6l/2) ) 



t_ / cos(0/2) -sin(0/2)\ 
"2 - l sin(0/2) cos(0/2) I 



where iïi = 7^(7^,0), fij = 7^^(7r, 0),iÏ2 = n{9,0) and iï^ = 7^(6',0). The rotation 7^(é', (/)) is defined by Eq. (U). 
The matrix a denotes the NOT operation. If not all control qubits are in the state |e), then the gate performs on 
the target qubit the unity operator 1 = r\i R^l R2 Ri- Finally, we may write the truth table of the multi-qubit 
(controlled)'ï-iï gate as follows 



l^no) ISlTlg + l) 

l*no)|e,„,+ ,> 
|*yeí)|5m,+ i) 



|*no)|em,+i) , 

\^ y,,) {cos 9 \g,n,^,) - e~'^^ sine \e^^^,)) , 



(7.23) 



where \^no) and {"^yes) are defined in Eq. ( |7.18| ). The multi- qubit controUed-R (CROT) gate (FIG.|Ï|) is performed 
on cold trapped ions by applying the evoluti on o perator (7.19) for the multi-qubit CNOT gates and the corresponding 
operator for the single-qubit rotations [Eq. ([7.3|)]. 
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FIG. 18: Scheme of a multi-qubit ( cont rollcd)'^-_R quantum gate. The gate acts on g -|- 1 qubits with q control qubits and the ríiq+ith 
qubit is target. R is defined by Eq. (7.22 ) in the basis {\g)mg^i , |e}m,^i } of the target qubit. Ri, Jï|, R2 and iíj s-^e also defined by Eq. 
(7.25) in the same basis. 



If the preparation of a particular class of quantum states does not require the introduction of a relative phase shift 
(j) between the basis states \g) and |e), then a reduced quantum lògic network is sufficient (FIG.p!9|). In particular, 
the rotation R = a R2cr R2 on the target qubit conditioncd by the state of control qubits can be realized according 
to the following truth table 



l^no) lffT?ig+i) 
|*no)|e™,+ i) 

|*yeí)|5m,+ i) 



|í'no)|em,+ i) , 



yesi 



costí \g. 



smW e. 



™5 + 



(7.24) 



^ J/es/ |Ct?x,+i/ 



\^y^s){sv[ie\g„,^^,) -f cos6l|e,„,_,i)) . 

The results for the multi-qubit controlled-iï gates are compatible with the scheme proposed in Ref. 0, where the 
decomposition of multi-qubit CNOT gates into the network of two-qubit CNOT gates has been presented as well. 
However, this decomposition may require many elementary operations in a particular realization of quantum lògic 
gates. It seems to be more appropriate for some practical implementations of quantum computing to implement 



directly multi-qubit CNOT gates (see Sec. IX ) 
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FIG. 19: Scheme of a reduced multi-qubit {contrpllcd)'ï-iï quantum gate. The gate acts on g + 1 qubits with q control qubits and the 
níq^lth qubit is target. R is defined by Eq. (7.24). R2 and iïj are defined by Eq. (7.25) in the basis {\g)mq^i, \^)mq^i} of the target 
qubit. 



VIII. QUANTUM LOGIC NETWORKS 



In this section we present quantum lògic networks as efïective tools for the synthesis of quantum coherent superposi- 
tions of intcrnal atòmic states. We provide two particular networks, where both of them apply to an arbitrary register 
of qubits. The networks consist of single-qubit rotations, multi-qubit controUe d-NOT and multi -qubit controUed-i? 
gates. Their implementation on cold trapped ions is described in detall in Sec. VII D and VII E . The generation of 
nonclassical motional states of a trapped ion experimentally is described in Ref. |73f~^ 

We keep the notation \g) and |e) for the logical states of the qubit also in this section. Firstly, let us introduce the 
network for the preparation of a totally symmetric state (with respect to the permutations) of N qubits, such that 
all qubits except one are in the excited state 



I*) = 



1 



\gee. 



\ege. 



\eeg. 



+ \eee ...g) 



(8.1) 



It has been shown that the maximal degree of bipartite entanglement measured in the concurrence is equal to 
2/N and is achieved when a system of N qubits is prepared just in the state (8.1). The synthesis of this state realizes 
the network in FIG.|^ assuming all qubits to be initially prepared in the state |e). The rotations Qj are given as 
foUows 



For more detalls we refer to our original paper 
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■2) 



2 
3 



Q2 



1 — r 
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FIG. 20: The network for the synthesis of the symmetric entangled state (í.l). The rotations Qj are given by Eq. (3.2). A'^ qubits are 
assumed to be initially prepared in the state |eee...e). 
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FIG. 21: An array of networks fo r th e synt hesis of an arbitrary pure quantum state ( p.3[ ) on three qubits. The initial state is Iq oq) and 
the rotations Uj are given by Eq. (8.4)-(B.6). The networks (a)— (g) generate gradually the respective terms in the superposition (3.3). 
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FIG. 22: A compact form of the array of the networks shown in FIG. pï[ 



Secondly, we propose an array of quantum lògic networks for the synthesis of an arbitrary pure quantum state for 
iUustration depicted on three qubits. However, the scheme is quite easily scalable on N qubits Let us assume a 
general state of three qubits in the form 



+ e''^^ai\gee) + e'^^^aslege) + e'^^^aelees) + e^'^'' a^\eee) 



(8.3) 



The state (8.3) can be realized by applying the array of the networks in FIG. 21 (shown in a more compact form in 
FIG.|2^) on the initial state \ggg), i.e. all three qubits in the state \g). We have denoted the rotations Uj as follows 



j = 0, . . . , 6 , 



(8.4) 



where aj = cosOj and bj = sin^j. The state (8.3) is given by 14 real parameters and the network in FIG. p2| preparing 
this state is also determined by 14 parameters (seven rotations), where 



and 



bn = a/1 



^0 ' 



A:=0 



j = l,...,6 



J = 1, . 



(8.5) 



(8.6) 



Thus, the mapping between the state under preparation (8_^) and the network (FIG. |22| ) is clearly defined 
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IX. SPEED OF QUANTUM GATES 



One of the most important requirements for thc implcmentation of quantum lògic (Sec. I) on a particular candidate 
quantum system is the physical realization of quantum gates on time scales which are much shorter compared to time 
scales associated with deco l·ier ence effects. We have represented quantum gates on cold trapped ions with unitary 



evolution operators ([7.3|), (7.5) and ( [7.q ) associated with làser pulses on the carrier and on the first red sideband. 
However, these operators are vàlid only in the Lamb-Dicke and weak coupling regime. Taking into account the 
mplete Hamiltonian ( 4.15| ) we have to deal with resonant and off-resonant transitions accompanied with Stark light 



co: 



shifts of the energy levels (Appendix ^). A detailed treatment of this problem was presented by Steane et al. in 
Ref. and we adopt some of their results in this section. The analysis of the speed of gate operations in ion traps 
was firstly discussed by Plenio and Knight in Ref. At first we discuss the speed of single-qubit rotations 

and two-qubit CNOT gates. Then we include some estimations for thc speed of multi-qubit CNOT gates with cold 
trapped ions. 



The single-qubit rotations are associated with the transition on the carrier (7.3) with the duration Ta — £tt/\X\, 
i.e. a ^TT-pulse on the carrier applie d on a given ion. However, we have to consider rather the evolution operator 
corresponding to the Hamiltonian ( 4.15| ) when we want to discuss unwanted off-resonant transitions. Directing 



the làser beam such that it is perpendicular to the z axis, the Lamb-Dicke parameter rj becomes equal to zero 
and off-resonant transitions |e)|r7,) ^ 1.9)1?^ ± |^|) for A; ^ do not appear in the dynamics. Therefore, one can 
make the làser coupling constant |AI large without the restriction on the weak coupling regime characterized by 
the condition |A| ^ v (Appendix H). We can assume |A|/27r ~ 300kHz, then we get typically for a 7r/2-pulse 
on the carrier Ta ~ 1 /is. 



The two-qubit CNOT gate (7.7) is realized by two 7r/2-pulses on the carrier [yím2^(0), yïma (t)], two 7r-pulses 



{Bl^ J and a single 27r-pulse {B„i J ) o n thc first red sideband. A pulse on thc first red sideband is represented 



with the unitary evolution operator ( 4.25 ) for k = —1. However, it was derivcd for an ideal case when off- 
resonant transitions and Stark light shifts were not considered. We can correct for the light shifts by tuning the 
làser on the frequency ujl — í^o — + Aw, where Auj corresponds to the light shifts caused by the presence 
of the carrier transitions. The imprecision caused by the excitation of off-resonant transitions can be corrected 
by applying a correction làser pulse with a correspondingly adjusted phase. This was accomplished in Ref. [T^ 
and the limit for thc duration of the operation corresponding to thc 7r-pulsc on thc first red sideband is given as 



where e = ^1 — F is the imprecision defined via the fidelity F, Er = S^K^/2m is the recoil energy of a single 
ion of the mass m, k= (27r/A) cosí5, A is the làser wavelength, h = 2nh and is the axial trapping frequency. 
The limit for the duration of the operation corresponding to the 27r-pulse on the first red sideband is the double 
of the expression given by Eq. ( |9.l| ). 

In TABLE || we give the estimations for Calcium ions ''°Ca^. We assume the angle between the làser beam and the 
z axis to be i5 = 60°, the làser wavelength is A = 729 nm and the axial trapping frequency is uJz/2tt ~ 700 kHz. Then 
we get the recoil frequency E^/h^ 2. 33 kH z and r/ = ^jEr/hujz — 0.06. 

The multi-qubit CNOT on Q ions ( 7.19| ) differs from the two-qubit CNOT gate only in the number of làser pulses 

required for its realization. The multi-qubit CNOT gate corresponds to two 7r/2-pulses on the carrier [vAm^^i(O), 
.4mg+i(7r)], a single 27r-pulsc on the first red sideband (;B^^U and {2Q — 2) 7r-pulses also on the first red sideband 
{Bl,[[, Bfnl^ , . . . , B^i^J). Thus, it requires all together 2Q + 1 làser pulses, where Q = q+1 rcfcrs to thc number of the 
ions involved in the gate (q control ions, one target ion). Then, in the spirit of the previous discussion the minimal 
total time for the realization of the multi-qubit CNOT gate on Q ions reads 

T = 2(T4 + QTb) , (9.2) 



where we assume T4 = 5 /iS and Tg is given by Eq. (9.f ) for a total number N {Q < N) of the ions confined in the 
trap. We give some estimations for the realization of multi-qubit CNOT gates in TABLE |[ where the number Q of 
the ions involved in the gate and the total number N of the ions in the trap are equal {N = Q). We stress this point 
because there is a difference if we realize a two-qubit CNOT gate [Q = 2) having just two ions in the trap {N = 2), 
then we get for the total time of the gate 

Ti = 2{Ta + 2T^=2) ~ 0.5 ms , (9.3) 
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iV 


Tb [l·is] 
F = 99% F = 75% 


T [ms] 
F = 99% F = 75% 


2 


124 24.8 


0.50 0.10 


3 


152 30.3 


0.91 0.18 


6 


214 42.9 


2.58 0.52 


9 


263 52.5 


4.74 0.98 


10 


277 55.4 


5.55 1.12 



TABLE I: A'^ is the total number of the ions confined in the trap and also the number of the ions involved in the realizati on o f the multi- 
qubit CNOT gate, Tg is the duration of the operation corresponding to the vr-pulsc on t he f irst red sideband given by Eq. ( |9.l| ) calculated 
for two vàlues of the fide lity ( F = 99%, F = 75%) and T is the total minimal time [Eq. (^.^)] for the realization of the multi-qubit CNOT 
gate on ions [see Eq. (r.l9)] evaluated for two different fidelities. 



or having a larger register of ten ions in the trap [N = 10), what gives 

T2 - 2{Ta + 2T^=^°) ~ 1.1 ms (9.4) 

at the same íidehty F = 99%. 

Any multi-qubit gate on a register of size N can be decomposed into a network of single-qubit rotations and two- 
qubit CNOT gates 0. However, there might be a possibiUty to reahze this multi-qubit directly, if a given physical 
system allows it. For instance, the multi-qubit CNOT gate on six qubits can be decomposed into the network of 12 
two-qubit CNOT gates including three additional auxiliary qubits 0|. In the case of cold trapped ions it requires the 
total time for the realization of the whole network {Q — 2, N — 9) 

Tg = 12 X 2(r4 + 2r^=^) ~ 12.7ms. (9.5) 



However, the direct implementation ( 7.l·E ) would take only {Q — N — 6) 

r4 = 2(T4 + 6r^=^) ~ 2.6ms, (9.6) 

which is about five-times less than the formcr case. We have again assumed almost the perfect fidelity F = 99% of the 
operation. We conclude that there can be quantum systems that may support a direct implementation of multi-qubit 
gates, rather than their decomposition into fundamental gates, what may bring advantages at experimental realization 
as well as in quantum state synthesis [[75[ . 



X. DISCUSSION 



A. Decoherence 



Throughout the paper we have discussed many aspects of cold trapped ions for quantum computing but we have not 
dealt with the decoherence which appears to be a main obstacle in achievements of experimental quantum computing. 
The reason was that our main goal has been togive a bàsic review and the discussion on decoherence sources and 
effects would refer more to an advanced study |7^, 77, Q Nevertheless, for the sake of completeness 

we would like to mention on this place some decoherence aspects met in the lab. We will foUow a detailed study of 



experimental issues in quantum manipulations with trapped ions given by Wineland et al. |22 . The decoherence will 
be met in a more general usage of this term. Thus, by the decoherence we mean any effect that limits the fidelity (the 
match between desired and achieved realization). Further, we will distinguish three categories. 

• Motional state decoherence is the most troublesome source of the decoherence in ion trap experiments and refers 
to the relaxation of two vibrational states |n = 0) and |n = 1) of a given motional mode used as the quantum 
data bus. The ions are cooled to the ground motional state \n — 0) and the excitation to the \n = 1) state is 
used for the transfer of Information on a distinct ion. However, this scenario is not ideal for several reasons: 



o Instability of trap parameters. We are simply not able to control all voltages as they undergo fluctuations 
and dephasing. 
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o We also have to count on (i) the micromotion, (ii) the Coulomb repulsion between the ions making the 
motional modes (except the COM mode) anharmonic in reality and (iii) stray elèctrode fields causing 
possible excitations of the ion motion. 

o We have considered just a single motional mode in our approach, but there are also other 3N — 1 modes 
present and the cross-coupling between the modes appears. If spectator 3A^ — 1 modes are not cooled to 
their ground motional states, the energy can be transferred to the mode of interest. This happens because 
the trapping potential is anharmonic in real and these higher anharmonic terms are responsible for the 
cross-coupling. 

o We should mention also inelastic and elàstic coUisions with the background gas, even though experiments 
are carried out in an excellent environment {p ~ 10^* Pa). 

• Internal state decoherence corresponds to the evolution when a pure state of the ion ji/í) — a\g) + f3\e) transforms 
into a mixture p = \a\'^\g) {g\ + \(3\'^\e) {e\. The ions demonstrate internal decoherence times of the order of seconds 
(Calcium) up to minutes and hours (Beryllium) . The type of the decoherence discussed here can be eliminated 
by a proper choice of metastable excited states with long lifetimes. 

• Operational decoherence refers to the precision of coherent laser-ion manipulat ions. There are several aspects 
that we have to consider: 

o When the ion is illuminated with a làser beam, one has to control the pulse duration and the phase 
adjustmcnt in order to avoid the preparation of unwanted states. 

o Duc to the làser spatial intensity profile, there is a probability (if the ions are spaced too closely) that the 
state of a neighbouring ion will be affected. 

o If we consider the standing-wave configuration we have to take care of the precise position of the ion in the 
node or the antinode of the standing wave what seems to be very troublesome. 

o Finally, off-resonant transitions are always present and we have to control the làser power very carefuUy to 
avoid their excitations. 

However, there is a way to eliminate the effect of the decoherence. We can encode Information into a decoherence-free 
subspace whose states are invariant under coupling to the environment Is^, Is^, Is^ . 



B. Ion trap systems 



We have been discussing cold trapped ions so far. Colà refers to the fact that all motional modes have to be cooled 
to their ground motional states because the dynamics assumes the precise control over the motional state. However, 
there have appeared other proposals referring to warm or hot trapped ions which assume an arbitrary motional state. 

• Poyatos et al. |Q proposed a scheme for the realization of two-qubit CNOT gates between two trapped ions 
using ideas from the atòmic interferometry. They split the wavepacket of the control ion into two directions 
depending on its internal state with a làser pulse. Then they address one of the wavepackets of the target ion 
changing conditionally its internal state and finally they bring together the wavepackets of the ions using another 
làser pulse. The ions communicate through the Coulomb repulsion and under ideal conditions the scheme is 
independent on the motional state of the ions. 

• Milburn et al. [s^ described two schemes for manipulations with warm trapped ions. Firstly, they use the 
adiabatic passage for the conditional phase shift, i.e. the phase of the ion is flipped if the motional mode is 
in the superposition of odd number states and the ion is excited. The COM mode in an arbitrary vibrational 
state is used for quantum lògic but all other motional modes are assumed to be cooled to their ground motional 
states. Secondly, they apply the idea of the coUective spin Q for faster gates. This idea has been also used for 
the introduction of multi-qubit gates for quantum computing [ p6| . 

• S0rensen and M0lmer |8^ proposed a novcl scheme bascd on the idea of bichromatic light (wi, tt'2)· Realizing 
the two-qubit CNOT gate they illuminate two ions with the bichromatic light coupling the states \gg)\n) and 
|ee)|n). They choose detunings far enough from the resonance with the first red and blue sidebands such that 
the intermediate states |eg)|n± 1) and |ge)|n± 1) are not populated in the process. The scheme is not sensitive 
on fluctuations of the number of phonons in the relevant motional mode. It is also possible to illuminate with 
the bichromatic light more ions and generate a multiparticle entangled state. Actually, these experiments were 
already realized in NIST and they generated the GHZ state with two ions \^p2) = 1/V2(|55) — ï|ee)) with a 
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fidelity F = 83% and also the GHZ state with four ions \tp4) = í/y/2{\gggg) + i\eeee)) with a fidelity F = 57% 
using Beryllium ions and the Raman scheme. Jonathan and Plenio proposed light shift induced quantum gates 
for trapped ions insensitive on phonon number in motional modes (thermal motion) ||8^ . 

• Finally, there has appeared a proposal of a scalablc quantum computer with the ions in an array of microtraps 
by Cirac and ZoUer |£Ï| detailed in Ref . ||9^ . The ions are placed in a 2D array of independent ion microtraps 
[ p3[ and there is another ion (head) that moves above this plane. If we position the head above a particular 
ion from the array and switch on the làser in the perpendicular direction, wc can realize a two-qubit gate. This 
operation allows us to swap the state of the ion to the head which can be moved immediately above a distinct 
ion in the array and transfer the Information onto it. The ions oscillating in the microtraps are not assumed to 
be cooled to their ground motional states. However, their motion can couple to the environment. It becomes 
relevant during the time when the ion interacts with the head but not in the casc when the head moves. 



XI. CONCLUSION 



In this paper we have tried to review achievements accomplished in the field of cold trapped ions. In the first part 
we have discussed in detail the ion loading and trapping process, the coUective vibrational motion of the ions. We 
have also given a detailed derivation of the Hamiltonian governing the dynamics of the system including the discussion 
of weak coupling and Lamb-Dicke regime. Further, we have reviewed làser cooling techniques and a detection process 
with experimental illustrations on Calcium ions. In the second part we have discussed the implementation of quantum 
computing using cold trapped ions. In particular, we have described how to realize single-qubit, two-qubit and multi- 
qubit quantum lògic gates. Finally, we have estimated the speed of quantum gates with cold trapped ions. The aim 
of this paper is to give an introduction to this field with many references on relevant papers and studies. 
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APPENDIX A: WEAK COUPLING REGIME 



The expression ( 4.15| ) corresponds to the complete Hamiltonian in the sense that it includes also off-resonant 
transitions. For instance, even though a sufíiciently intense làser is tuned on the carricr, the off-resonant transitions 
on the sidebands are also present and they cause imprecisions and pcrturbations in the dynamics. This can be avoided 
by setting the làser intensity / oc |Ap, i.e. the làser coupling constant [Eq. ( 4.12| ) or ( 4.13| )], sufíiciently small. In what 
foUows we will determine conditions (characterizing the weak coupling regime) under which we can neglect off-resonant 
transitions. 



1. Off-resonant transitions 



The implementation of quantum gates on cold trapped ions requires làser pulses on the carricr [FIG.^(a)] and on 
the first red sideband [FIG.|(b)]. Therefore, we will discuss the dynamics on these two spectral lines. At first, let us 
assume that the làser is tuned on the carrier ((5 = 0). The closest off-resonant transitions (detuned by the frequency i^) 
are on the fir st bl ue and on the first red sideband [FIG.|2^(a)]. We will consider only the respective terms in the 
Hamiltonian (4.15) and drop down the index j. In the Lamb-Dicke limit we can write 



Hi = —a+ 



hx 



hX 



-{iri)(j+à e 



H.c. 



(Al) 



Further, we assume the ion to be initially in the state |í7)|?t-), we apply the Hamiltonian (Al) and calculate the 



probability of off-resonant transitions. The dynamics governed by a time-dependent Hamiltonian 7í(t) is described to 
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the first order by the unitary evolution operator 

Ü{t,to) 



n{t')df . 



(A2) 



to 



Then the probability to find the ion (initially prepared in the state |g)|rí,)) in the state \e)\n + 1) (i.e. undergoing the 
ofï-resonant transition on the first blue sideband) is 



Pr 



(e|(n + l|í7i|.g)|n)|^ 



\X\^rj\n+l) 



■ sm 



iy{t - ío) 



(A3) 



where Ui is given by Eq. (A2) for the Hamiltonian (Al). The probabihty to find the ion in the state \e)\n — 1 
corresponding to the ofF-resonant transition on the first red sideband is given as 



Pn^\{e\{n-l\U^\g)\n)\^ 



v{t - to) 



(A4) 



If thcrc is no population transferred via the ofF-resonant transitions to the states |e)|n+ 1) and |e)|rí — 1), i.e. Pb <C 1 
and Pfí <C 1 at any time í, we can neglect these ofï-resonant transitions. Then we get the conditions of the weak 
couphng regime for the transition on the carrier in the form 

(A5) 
(A6) 



and 



\X\T]\/n+ 1 < 1^ 



We can also avoid the off-resonant transi tion s by sctting the làser beam perpendicular to the z axis {■& = 7r/2). 
Then the Lamb-Dicke parameter [see Eq. (4.9)] is equal to zero (rj ~ KCOS'dzo) and the coupling on the off-resonant 
transitions vanishes. 



5 = 



le>ln> 



le>ln-l> 



T|A, 





X 







lg>ln-l> 



lg>ln> 

(a) 



le>ln+l> 

— y 

lg>ln+l> 



5 = -v 

le>ln-l> 



le>ln> 
1— 



r|À. 



lg>ln-2> 



lg>ln-l> 

(b) 



lg>ln> 



FIG. 23: (a) The làser is tuned on the carrier (5 = 0) with the coupling constant A. The off-resonant transitions are present on the first 
blue and red sideband with the coupling constant rjX. (b) The làser is tuned on the first red sideband (5 = —u) with the coupling constant 
r}X. The ofï-resonant transitions are present on the carrier with the coupling constant A and on the weak second red sideband (rj·^A). 



Analogically, we can assume the làser to be tuned on the first red sideband {6 — —v) and the closest off-resonant 
transitions are on the carrier and on the second red sideband [FIG.|2^(b)]. However, the strength of the second red 
sideband is of the order of jf' and we can omit it in the Lamb-Dicke limit. Then the respective Hamiltonian is given 
as 



hX 



-{iri)à+à 



hX 



H.c. 



(A7) 



The probability to find the ion (initially in the state in the state |e)|n) (after the ofF-resonant transition on the 

carrier) can be calculated as 



Pc = |(e|(n|í72|.g)|n)p = ^sin 



iy{t-to) 



(A8) 



where C/2 is given by Eq. (A2) for the Hamiltonian (A7). We can neglect the off-resonant dynamics if Pc ^ 1 at any 
time t. Then for the weak coupling regime on the first red sideband applies 



|A| < I/. 



(A9) 



Even though the transition on the carrier is ofF-resonant, it has stronger coupling in the Lamb-Dicke limit than the 
first red sideband. Therefore, it is very important to foUow in the experiment the constraint given by Eq. (A9). 
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2. Stark light shifts 



Besides the population of thc ofF-resonant levels there is another source of imprecisions in the state manipulation. 
However, it is weaker and it doesn't require any special constraints on physical parameters except those for the weak 
coupling regime. When the làser drives a transition between two levels there appears a frequeney shift called Stark 
light shift caus ed by the presence of other spectator levels. Therefore, in the experiment we have to consider the 
dctmiing ( 4.14| ) rather in the form 



S = LÜL — LÜQ — Alü . 



(AIO) 



where Au corresponds to the Stark light shift. Thc highcr thc làser intensity, thc more significant the light shift is. 
We can correct for this effect by shifting the làser frequeney (wl + Auj) which tunes thc transition back to the 

resonance. Further, we will estimate the frequeney shift Alü. 



Lct us consider the Hamiltonian ( |4.4| ) and transform it to the interaction picture ( Í.IC ) with Uo = exp{—iHot/h), 
where Hq = {hujL/2)àz· Then we get 



Tí = Tío + V = (^-yà, + hi^à^à^ + (^ya+ e^''("+^') + H.c.) , 



where in the Lamb-Dickc limit the interaction tcrni reduces to 



(All) 



(A12) 



with the first term corresponding to thc transition on the carrier, the second term to thc first bluc sideband and the 
last one to the transition on thc first red sideband. In the second order of the time-indcpendcnt perturbation theory 
(the first order gives no contribution) we can write for thc shift of the energy levels 



AE, 



19)1' 



E 



(ff|(n|V|e)|m)(e|(m|V|.g)|n) 
h6 + hv{·n — m) 



(A13) 



and 



AE, 



\e)\r 



E 



(e|(n|V|ff)|m)(g|(m|V|e)|n) 
—hS + hv{n — m) 



(A14) 



For instance, for the transition on the first red sideband {5 = —v) we can calculate 



AE\ 



|9>I»> 



Av 



AE, 



\e)\n-l) 



and the corresponding light shift can be estimatcd as 

= (^Aí;|e)|„_i) - AS|g)|„)^/ri; 



Aiy 



2v 



(A15) 



(A16) 



If we choosc |A|/27r = 50kHz and i//27r = 700 kHz so that the condition ([A9|) holds, then Auj/2tt ~ 1.8kHz. 



APPENDIX B: LAMB-DICKE REGIME 



In the relation for thc coupling consta nt ( 4jl) we can expand thc exponential function to thc Taylor series about 
the value rjj — and usc thc expression ( 4.22 ) for the Laguerre polynomial. Then we get 



\k\ 



n+\k\ 



(n+ |fc|)! 

71 H 

n 



(Bl) 
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where 0{i]'^) denotes the terms proportional to the fourth and higher powers of rjj. In the Lamb-Dicke regime we 



consider the dependence of ííj' on the parameter rjj only to its lowest order, i.e. in Eq. (Bl) we neglect all terms with 

any higher power than r/^'^^ and we get the coupling constant Í2"''^ given by the expression (4.29). This approximation 
can be done only if the conditions 



2 

< 1 (B2) 



and 



2 



ín+\k\\ ín+\k\\ , n , , 



are satisfied. We will refer to the condition (|fc| = 0) 



Vj^{n) + l<t:l (B4) 

as the Lamb-Dicke limit, where (n) is the average number of phonons in the respective vibrational mode . 

• The Lamb-Dicke limit corresponds physically to the situation where the spatial extent of the vibrational motion 
of the ion zq is much smaller than the wavelength A of the làser, where rjj ~ kzq and k = 27r/A [see def. in 



Eq. (4.9)] 



• The Lamb-Dicke limit can be physically interpreted also from a different point of view. We may rewrite the 
Lamb-Dicke parameter of the j ion of N ions in the COM mode to the form rj^ = Er/hujz, where Er = fi^K^ /2mN 
is the recoil energy. It can be shown that the trapped ion emits spontaneously photons of the average energy 
Hujo — Er, where luq is for the atòmic frequency. Taking into account the Lamb-Dicke limit (i?,. <^ hcuz) we may 
say that during the spontaneous emission the change in the vibrational state of the ion is very unlikely. In other 
words, the trapped ion in the Lamb-Dicke regime decays spontaneously mostly on the carrier (wq ~ luq — Er/h). 



The Lamb-Dicke parameter for a single trapped ion equals to r] and for an ion from the string of N ions in the 



COM mode is given as rjj = rj/yN. It means that we can reach the Lamb-Dicke limit (B4) for N ions even if 
the limit is not fulfilled for single ions 
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